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Thickness of clogging layer at left boundary , 
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c Total resistance to inter-aquifer flow , 

c Resistance to int er -aqul f er flow due to unconfined aquifer ; 
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d Average depth of flow In unconfined aquifer , 

e Potential rate of evaporation , 
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Resistance to moisture flow lii soil , 
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Time from begining of first rainfall to the time when 
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Duration of first rainfall , 

Time from begining of first rainfall till soil surface 
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second rainfall starts , 

Nondimensional time ; 

Actual rate of transpiration , 

Potential rate of transpiration ; 
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Potential rate of transpiration per unit depth of root 
zone , 

Recharge from top , 
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U Nondimenelonal rechaige from top , 

X Horizontal distance measured from origin , 

X Non-dimensional horizontal distance ; 

y Vertlcle distance measured from the datum , 

z Vertlcle distance measured from ground surface , 

Z . Non-dimensional vertical distance measured from ground 
surface , 

Index used in defining transpiration rate - effective 
saturation relationship , 

5 Diffuaivlty of unsaturated soil , 

r/ Effective porosity , 

Effective saturation , 

* 

Effective saturation corresponding to threshold moisture 
content , 

Effective saturation corresponding to wilting point , 
wp 

y Specific weight of soil water , 

6 . Volumetric moisture content , 

Q Volumetric moisture content at saturation ; 

s ’ 

K Pore size distribution index , 

p Conductance of the semi-pervious layer ; 

‘ A macroscopic physical property of soil which indicates 
the amount of work doite per unit volume of soil reiiuited 
to drain a saturated soil to the wilting point ^ 
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CHAPTER 1 


INTRODUCTION 


1.1 INTRODUCTION 

Increased demand of water supply for a variety of purposes 
require conjuctive use of surface and ground water resources 
Recently lot of emphasis has been put on ground water resources 
since the surface water resorces are limited and less depend- 
able Ground water is available in a number of water bearing 
formations beneath the surface of earth called aquifers The 
topmost aquifer usually occurs with a pheratlc surface and is 
called an unconfined aquifier Flow in unconfined aquifers is 
important not only from the point of view of extraction of 
water but also from agricultural point of view In fact it is 
this aquifer which is subjected to maximum f luctutations Some 
of the problems involving uneonfined flow are 

(i) Flow between a river and an unconfined aquifer with or 
without recharge to the aquifer , 

(li) Flow through and below an earth dam , 

(ill) Flow to a well in an unconfined aquifer ; 

(iv) Flow between two parallel drains or ditches provided for 
controlling ground water , 

(v) Flow in an unconfined aquifer overlying a semi-pervious 
strata ; 

(vi) Artificial recharge problems 

All the problems involve finding of water table heights as a 
function of apace and time. Though man la concerned with uncon- 
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fined aquifers for a long time analytical solutions for only a 
few problems exist even to date in contrast with the confined 
aquifers This is mainly because of the following reasons 

(a) The partial differential equation governing the ground water 
flow in unconflned aquifer is non-linear 

(b) The boundary at the top of unconflned aquifer is not known 
until a solution is obtaitied , and the solution depends on the 
boundary whereas the flow domain in a confined aquifer is 
initially known 

The exact analytical solutions for many complicated ground 
water flow problems are still not available owing to the above 
mentioned difficulties Numerical methods are usually employed 
in such cases. Analytical solutions, if available are always 
preffered, even if they are often based on simplifying assump- 
tions Analytical solutions provide following advantages 

fa) It is convenient to use analytical solution, if one exists, 
(b) They are useful in verifying numerical methods 
fc) The effect of various parameters can be studied in a quick 
and easy way with analytical solutions 
(d) Also such solutions are useful In priliminary studies of 
aquifer developement schemes 

Various approaches have been used by investigators to find 
analytical solutions Theory of equipot ent ial flow has been used 
extensively for finding analytical solution of two dimensional 
problems. Conformal transformations are used to simplify the 
flow situation But this approach has serious limitations for 
analysing transient flow situations 
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The approximate analysis of shallow flow in unconfined aqui- 
fer over an impermeable bed was initiated by Dupuit by maklne 
the well known Dupuit 'a assumptions These asaumptiojis t educe 
2-D flow situation to 1-D and 3-D flow situation to 2-D 

Eoussinesq derived the governing equation for the unsteady 
state flow over sloping fully impervious bed with the help of 
Dupuit ‘s assumptions Additional terms are required to be added 
to the original Boussinesq equation to take into account rain- 
fall recharge and recharge from or leakage to the underlying 
aqul f er 

To overcome the difficulty imposed by nonlinearity, various 
types of linearisation have been proposed in the literature. 
The linearised Boussinesq equation provides a powerful tool for 
treating flow problems in unconflned aquifer in order to find 
an engineering solution Once the equations are linearsed the 
principle of superposition can be utilised Such linearisation 
is of much value for treating transient flow problems. 


When the underlying strata is not fully Impervious, leakage 
from or to the aquifer takes place In addition to recharge from 
top, if present. The magnitude of the inter-aquifer flow will 
depend upon the difference of heads in unconfined aquifer and 
underlying aquifer, and in general this will be nonuniform. The 
recharge from top however is usually taken to be uniform If the 
the piezometric surface In the underlying aquifer is too high, 
artesian condition may occur Usually such condition occurs in 
the following circumstances . 

(a) In vally bottomlands where alluvial fans are often under- 
lain by lensed semi-pervious formations 
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In coastal plains where the high oceti level causes auch 


(b) 

problem Typical flow situation at a Dutch polder la of 
this type 

(c) At foundation construction sites 

(d) Uhen land Is protected by river levees and the artesian 
pleasure la Imposed by navlgatloii Improvements or periodi- 
cally imposed by flood levels 

1.3 OBJECTIVE 

The purpose of this study Is to Investigate flow problems In 
unconflned aquifer on sloping semi -pervious beds Following 
types of problems have been considered in this context: 

(1) An exact solution for the steady state flow over horizontal 
semlpervious bed when artesian head line is also horizontal 

(2) An exact solution for the unsteady state flow over horizon- 
tal barrier between parallel drains provided on the barrier 

(3) A linearised solution for steady state flow over sloping 
bed when artesian head line Is also parallel to the bed. 

(4) A general solution for unsteady state flow for 


(1) 

Trans! ent 

r echarge 

from top 


Cil) 

Leakage from or to 

unconflned 

aqui f er 

(ill) 

Arbitrary 

initial c 

ondition 


(iv) 

Di f f er ent 

boundary 

conditions 


(V) 

Horizontal 

as well 

as sloping 

bed 


1.3 APPROACH 

Methods of Laplace Transform and Fourier Transform have been 
extensively used for solving partial differential equations. 
Many investigators have used this approach to find valuable 
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solutions for ground water flow problems. But in the present 
study, method of substitution la used Simplified equation thus 
obtained can be easily solved for a wide variety of boundary 
conditions In case of transient flow problems the substitutions 
will make the boundary conditions transient In these problems 
Duhamel's theorem will be used to deduce solution for transient 
boundary conditions. The transformed initial condition can be 
easily Incorporated using Fourier series Then a process of back 
substitution will lead to the desired results 

The above mentioned approach has the advantage that only a 
simple partial differential equation requires to be solved for 
which solution exists in literature. Again this approach la also 
applicable to problems of other fields where similar partial 
differential equation is required to be solved e.g. solute 
transport problems 
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CHAPTER 2 


LITERATJRE REVIEW 

Literature connected with flow towards drains and streams ia 
reviewed in a chronological order in the following section : 

The analysis of flow in unconfined aquifer was initiated by 
Dupuit (1863) by making the assumptions that the flow lines are 
horizontal and hydraulic gradient at any section ia equal to 
the slope of pheratic surface at any section These assumptions 
reduce a 2-D problem to 1-D and 3-D problem to 2~D 

Using Dupult's assumptions Boussinesq (1904) derived the 
partial differential equation governing unconfined ground water 
flow in non-def ormabl e homogeneous porous media lying over a 
sloping impervious bed 

Hooghoudt (1940) used Dupult-For chhelmer approach to find an 
exact solution for steady state drainage of uniform rainfall 
by parallel drains in a horizontal aquifer To account for con- 
vergence he proposed the use of equivalent depth instead of 
actual height of drains above Impervious barrier 

Muskat (1946) found a solution for flow between two parallel 
drains over horizontal semi-pervious bed by conformal mapping 
He assumed that the horizontal line passing through the maximum 
height of the water table is a stream line^ drain radius was 
negligibly small and flow near drains is essentially radial. 
Recharge from top or the presence of any seml-pervloua layer at 
bottom was not considered. 
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Uerner ( 1953, 1957) Introduced rechar£e term In the 

2 

Bouaalneaq equation and uaing tranaf ormation h =y, linearlaed 
the equation and aolved it £or flow between parallel reaervolra 
and for flow between ground water divide and a reaervoir ualng 
Laplace tranaformation approach 

Glover (1953,1954,1960) linearlaed the Bouaalneaq equation 
uelng average depth of flow and found analytical aolutlon for 
the caae of flow between parallel dralna over impervioua beda 

Klrkham (1958) conaldered the problem of aeepage of ateady 
rainfall through aoil into dralna aa a boundary value problem 
and applied the potential theory approach to find a aolution. 

Kraljenhoff van de leur (1958,1962) aolved the Bouaalneaq 
equation for predicting ground water table between two parallel 
dralna in a horizontal aquifer aubjected to conatant recharge. 

Waaaland (1959) waa one of the firat to treat the problem of 
tranaient recharge analytically He approximated the input terma 
by a plecewlae conatant function and uaed the aolution for 
conatant rate of recharge to conairuct the aolution for tran- 
sient recharge 

Schmid and Luthin (1964) found implicit aolution for ateady 
state flow between parallel drains over sloping barrier uaing 
method of subatitutiona Spacing of drains was found explicitly 
The maximum height was compared with the experimental aolution 
obtained with the help of Hele-Shaw model and was found to be 
sufficiently accurate 

Moody (1966) uaed numerical approach to provide correct 
solution in graphical form for flow between two parallel dralna 
over a horizontal impervious barrier The eonvergece to drains 
was accounted for by uaing Hooghoudt’e equavalent depth for 
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which he developed ftimpler expree^ilone The error due to uae of 
simpler expressions was found to be leas than 

Hlnesly and Kirkham (1966) considered the flow between two 
parallel drains over horizoi tal semi-pervious bed subjected to 
a steady recharge as a boundary value problem which they solved 
with the help of Fourier Transform The head loss of water flow- 
ing above the drain level was neglected and drains were assumed 
to run half full The maximum height of water table was taken 
as the head midway between the drains 

Luthin and Guit Jens (1967) studied the transient water table 
falling between parallel drains in a sloping aquifer with the 
help of a Hele-Shaw model They gave empirical relations to 
estimate drain spacing and indicated that rate of fall of water 
table was almost independent of slope for slopes upto 30^ 

Chauhan et al . (1968) obtained an analytical solution for 
time varying end conditions for flow in a sloping aquifer. They 
also obtained experimental solution wth a Hele-Shaw model and 
found that the Boussinesq equation is capable of characterising 
transient water table for drains penetrating to the full depth 
Of impermeable layer upto 8% slopes 

E C Childs (1971) pointed out that for flow over sloping 
beds largely governed by bed slopes, streamlines will be para- 
llel to bed Instead of being horizontal , the equlpot entlal line 
being normal to it. Bach equipot entlal line is still labelled by 
the height of point where it cuts the pheratic surface measured 
from a horizontal datum Thus he obtained a new expression fbr 
flow using Darcy's law and found an exact analytical solution 
for flow between two parallel drains transverse to the bedslope 
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for the case of no recharge 

Marino (1973) considered a stream aquifer system over an 
impermeable base and derived analytical solutions for the cases 
when bank of stream was clogged due to the presence of a thin 
layer for both a emi- inf init e and finite aquifers The initial 
water table was taken to be horizontal and in case of finite 
aquifer the other boundary was taken to be of no flux type No 
recharge from top was considered and Laplace Transform was used 
to obtain a solution He (1974 a,b,c) also presented analytical 
solutions for transient pheratic surface in seml-infinte and 
finite horizontal unconfined aquifers which recleves localised 
recharge at a constant rate and discharges into a neighbouring 
reservoir The Initial water table was taken to be horizontal 
and Laplace Transformation was applied to find a solution 

G.D Towner (1975) employed proper substitutions in Childs’ 
(1971) differential equation for steady state flow between two 
parallel ditches over sloping bed with uniform rainfal 1 , thereby 
obtaining similar integrations as Schmid and Luthin Thus impli- 
cit analytical results were found for flows controlled largely 
by sloping bed 

Singh and Jacob (1977) used the method of functional trans- 
formation for linearisation of Boussinesq equation. The equation 
so linearised was solved for flow between parallel drains or 
streams with unsteady recharge for arbitrary initial conditions 
and transient boundary condition of Dirichlet type. Excellent 
agreement was found between linearised and numerical solutions 

Mualem and Bear (1978) proposed a new transformation and 
linearisation technique for steady pheratic flow over a sloping 
semi-pevious bed. The solution of the linearised equation was 
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verified with a Hele-Shaw analog Excellent agreement wae found 
between analytical and experimental reaulta 

T G Chapman (1980) gave a aimpler expression for discharge 
over a sloping impervious bed utilizing a new approximation iti 
Chi Ids ’ ( 1 9 71 ) expression The new approximation was found valid 
upto a slope of SO'’ 

Sakkaa et al (1981) fitted simple equations for Hooghoudt 
equavalent depth by using method of least squares and indicated 
that maximum errors were less than 1% 

Ral and Singh (1981) obtained a solution for transient water 
table in a horizintal seml-inf Init e aquifer subjected to loca- 
lised transient recharge linearly decreasing with time and then 
remaining constant using Laplace Transformation approach. 

Bruggeman (1983) considered the steady flow between parallel 
drains over seml-pervlous bed as a boundary value problem as 
Kinsley and Kirkhatn (1966) did, but some extra parameters such 
as the resistance of a thin layer at bottom, resistance against 
entry of water into drains and anisotropy were taken care of 
Uesseling and Uesseling (1984) used the Dupuit-Forehhemer 
approach to find the solution for flow between parallel drains 
over horizontal semi-pervious bed subjected to steady recharge 
The resistance of a thin loamy layer at bottom was taken care 
of Using Glover's type of linearisation they found a solution 
for steady as well as unsteady flow 

Tolikas et al . (1984) presented an approximate analytical 

solution for the problem of recharge of a semi-infinite aquifer 
from a stream with uniform initial condition and a step incre- 
ase of piezometric head on the boundary. The 1-D problem was 
reduced to an ordinary differential equation through Boltzman 
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Transformation and a technique exploiting some basic character- 
istics of the exact solution leads to an approximate polynomial 
solution of the problem 

Yates and Uarrick (1985) developed an exact solution for 
steady state flow in a sloping Dupuit aquifer subjected to a 
constant uniform recharge In case of leaky substratum they 
suggest that aritesian recharge should be added to the recharge 
from top assuming it to be uniform 

Latinopoulos (1986) developed an analytical solution for 
pheratlc surface in finite horizontal aquifers when one of the 
boundary condition was of third type, the other being either of 
Dirlchlet type or Newmann type. The initial condition was taken 
as horizontal water table and a periodic recharge ( constant or 
zero ) was considered to take place The solution was obtained 
using Laplace and Finite Fourier Transforms. 

Bazaraa et al . (1986) studied the aritesian and anisotripic 

effects on drain spacing. They indicated that large plezometrlc 
heads may necesslate small spacings which may be uneconomical 
The anisotropic effects were studied by converting analotropic 
system into a fietlous isotropic system keeping the vertical 
dimensions unchanged It was found that drain spacing Increases 
as anisotropy ( ) increases for impermeable beds For 

permeable beds and shallow depth of aquifers spacing increases 
Initially with anisotropy till magnitude of aritesian flux is 
reduced a level comparable to rainfall recharge. Thenafter it 
decreases . 

Uolters (1986) compared five different formulae given by 
Muskat (1946), Young ( 1 986 ) , Hlnesly ( 19 66 ) , Bruggeman (1983) and 
Uesseling and Uesseling (1984) for certain cases He indicated 
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that Ueaseling and Uessellng (1984) formula in combination with 
Hooghoudt equavalent depth gives equivalent results as Hinesly 
and Brugeeman, and maximum height of water table so obtained is 
higher than those given by Iluskat’s and Young's formulae 

Youngs (1986) used the concept of potential flow theory and 
conformal mapping to find out water table heights between two 
partially penetrating drains draining the artesian water in a 
horizontal aquifer The solution involves integrals which can be 
evaluated only by numerical methods. To avoid this difficulty he 
gave an empirical equation which fitted the correct solution 
well 

Ram and Chauhan (1987) found analytical solutions for rise of 
water table between two parallel drains in response to steady, 
linearly increasing and exponentially decreasing recharges with 
zero initial condition The analytical solutions were verified 
experimentally with a Hele-Shaw model for slopes varying from 
0-0 078 % Large errors were found near ends due to convergence 
of flow The analytical solution was found acceptable though the 
need of proper selection of average depth of flow in this case 
( when drains are on the barrier ) was indicated. 
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CHAPTER 3 


FORMULATION OF THE PROBLEM 


3.1 INTRODUCTION 

Generally ground water flow ptroblenia in unconfined aqulfera 
involve the determination of phreatic surface and its variation 
with respect to space and time Sometimes the rate of leakage to 
drains and streams Is also desired Some of the common flow 
situations are flow between two parallel channels la), 

flow between a no flux boundary and a river (/ftf-S.lb) and flow 
between two parellel pipe drains (/ftf-3 Ic) employed to control 
the water table The bed of the aquifer may be sloping or aemi- 
pervious or both 3 . Id , e ) . In case aquifer bed Is sloping, 

the phreatic surface will no longer be symmetrical and hence 
location of ground water divide cannot be predetermined as in 
the case of aquifer with horizontal bed In all the above cases 
a transient recharge may be taking place from top In practice 
such a recharge arises due to rainfall or due to excess irriga- 
tion or artificially applied water at the surface In addition 
to the recharge from top, flow of water to (/f^-S.ld) or from 
le) the unconflned aquifer may take place if the aquifer 
bed is not fully imper'^ious -Such situations are quite common in 
practice In a typical situation C/i-8“3 ld,e) most frequently 
encountered in practice, the unconfined aquifer is underlain by 
thin semi-pervious stratum which in turn overlies an aquifer If 
the piezometric surface in this semi-confined aquifer is higher 
than phreatic surface, recharge to unconfined aquifer will take 
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Pig. 3.1 (c) 

Fig. 3.1 Some ground water flow altuatlone 
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Plfezometrlc head line 
for underlylnfi aquifer 




Rainfall 


Phreatic surface 



Sloping aeml-pervious bed 



Underlying aquifer 



Fig. 3.1 (contd ) Some ground water flow situations 
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Rainfall 



Underlying Aquifer 



Fig. 3 1 (e) 


Fig 3.1 (dontd.) Some ground water flow aituationa 
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place This recharge may be appreaciable necealatlng provision 
of drains to eliminate this unwanted water If the piezometric 
surface is low, the unconfined aquifer will be losing water The 
recharge to or from unconfined aquifer in magnitude will depend 
upon the difference between the elevations of piezometric and 
phreatic surface and since either may not be uniform, it will be 
nonuniform Therefore the practice of adding it to the recharge 
from top, which is usually uniform, will lead to errors Thus 
recharge from top and bottom need separate treatment Special 
treatment Is also needed for streams with clogged bed bounding 
the aquifer Usually beds and banks of streams are overlain by 
a thin s emi -pervious layers called clogging layers (/lg“3.1e) 
Presence of such a layer will affect the flow rate to or from 
stream consldeiably and hence the water table heights 
3.2 GOVERNING PARTIAL DIFFERENTIAL EQUATION FOR FLOW IN 
UNCONFINED AQUIFERS 

The governing partial differential equation for ground water 
flow in an unconfined aquifer consisting of homogepeous and 
isotropic nondeformable porous media 
combining continuity equation with 
control volume shown in the 

dh (h -h) 

H O 

+ T) 

Sy. c 

Uher e 

2 

q flow rate in positive x-direction, L T , 

f} effetive poros ity. dimens ionl ess , 

X horizontal distance measured from origin, L; 
t time.T; 


can be easily derived by 
Darcy's law Consider the 
From continuity, 

« . C3.1) 
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Fig. (3. 2) Definition sketch for sstureted ground 
v&ter flow 


Id 




c 

h 

o 


h 

From 


realatance o£ underlying aemi -pervi ous layer, T, 

vertical distance to piezometric head in the underlying 

aquirer measured from bed,L, 

vertical height to phreatic surface measured from bed,L; 
tranaiant recharge from top,LT''\ 

Darcy's law 


q = -k h 

K m 


y 


4>X ' 


C3.2) 


where 


k saturated permeabi 1 i ty.LT" ; 

a 

h verticle distance to phr*atlc surface 

y 

datum,L ; 

I slope of bed , dlmenalonl ess 
From C/ig-3 2) 

h = h + y 

y 


measured from 


(3 3) 


Therefore , 


q = - k h (- 

H • 




dx 


+ I) 


. .( 34 ) 


Combining eqn (3 1) and eqn (3 4) we get the final equation as 


Sh 

n 

at 




ah 

h ( + 

ax. 



+ w 

c 


(3 5) 


3.2 1 ASSUMPTIONS INVOLVED 

The above derivation Involves following assumptions 

(1) The acquifer is homogeneous and isotropic. 

(2) The acquifer consists of nondef ormabl e media. 

(3) Darcy’s law is valid. 

(4) Dupl t-Forchheimec assumptions are valid 

(5) Recharge rate from top is small enough so that it gets 
refracted in the direction of groundwater flow itself. 
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(6) Flow of water In the underlying semipervl one , . le 
vertical and as soon as it enters the unconfined acquifer 
it gets refracted in the flow direction prevallnig in the 
unconflned acquifer 

(7) Bed slope is small 

(8) The flow boundaries are vertical 


3.3 INITIAL AND BOUNDARY CONDITIONS 

The solution of any partial differential eqation requires 
additional information in the form of boundary conditions to 
evaluate the arbitrary constants appearing therein For unsteady 
flow problems initial condition is also required to determine 
the solution uniquely 
3.3 1 INITIAL CONDITION 

Initial condition in its most general form can be written as 


h = f(x) 


at t = 0 


( 3.0 


If the field data is available, function f(x) can be evaluated 
by fitting a polynomial for all the points or by representing 
it by a series of parabolas each fitting three consecutive 
points of fCx) Thus in general the Initial condition can be 


taken as, 


f(x) . 2 ( (X) - '5 


(a^ X + b- X + c') 


C3.7) 




1 


x^ < X < x^^2 
1 = 1,3,5 ... 


Xj^ « 0 

^nt2 = ^ 


.. (3.8) 

. . , ( 39 ) 


In the solution initial condition is needed only to carry 
out an integration of the form 
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f(x)CA sin a X + B cos ot x ) dx 
332 BOUNDARY CONDITIONS 

When the water surface elevation above the bed la known 

at the boundary l,e in the neighbouring stream , the boundary 
condition is said to be of Dlrichlet type and ia given by 

h=h^ atx=0 (310) 

h=h atx=l .(311) 

s X 

If fully impervious bed ia preaent at the boundary the boundary 
condition is of no flux Neumann type and ia given by 


=0 atx=0 (3.12) 

dx 

If banka of the streams forming the boundary are clogged due 
to the preaence of a clayey material of low permeability and 
practically no storage, the boundary condition will be of Cauchy 
type and will be given aa 



+ k 

8 




k 

^ T- 


0 


0 


at X = 0 


at X = 1 

X 


(3 14) 

(3 15) 


All the three types of boundary conditions mentioned above can 
be represented as 


21 



elt X 


0 


C3 16) 


“k 


3 


b 

1 


dh 


dx 


b 

4 


+ k h = k h 

4 4 4 


k 


s 


b 

B 


dh 


b 


fi 


+ k h = kh atx=l 

BBS X 


(3 17) 


where 

thickness of clofifilng layer at left boundary, L; 
b^ thickness of cloegine layer at right boundary, L; 

permeability of clogging layer at left boundary , LT"*^ , 
permeability of clogging layer at right boundary , LT"* ; 
transient water surface in left channel, L; 

h transient water surface in right channel, L: 
s 

In ecin (3 16) = 0 gives rise to Dirlchlet type boundary 

condition whereas k^= 0 will give rise to no flux Neumann type 

bondary condition. In case none of them are zero it represents 
Cauchy type boundayr condition.Similarily , with b^=0 and k^a 0 
in e<;t*^-(3 17) give rise to Dirlchlet and no flux Neumann type 
boundary condition respectively. 

3.4 RESISTANCE AMD CONDUCTANCE 

The inter-aquifer flow between unconflned and semi - con f Ined 
aquifer has to experience certain resistance 'c' Reciprocal of 
this is termed as conductance ’p’ and is given by 

p = 1/c (3.18) 
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Thicker the aeparaling layer between two aquifers or lower Itfl 
permeability larger will be the vertical resistance to inter^ 
aquifer flow In absence o£ this layer the flow has to still 
encounter certain resistance due to the unconflned aquifer Thus 
the total realatance to the inter -aqui f er^ comprise of tvo parte 

c = c' + c . (3 19) 

u 

each part being given by 

c'= b'/ k' (3 20) 

c = d / k (3 21) 

u a s 

where 

b thickneee o£ the eemi-pervioue layer separating the two 
aqul f era , L , 

k permeability of the separating layer, LT”*, 

d average depth of flow.L: 

a 

k saturated permeability of unconfined aquifer, LT 
& 

c total resistance to inter-aquifer flow.T; 

c' resistance to Int er -aqul f er flow due to uneonflned 
aqui f er , T , 

c^ resistance to Inter-aquifer flow due to separating 
layer , T 

The vertical resistance resistance due to uneonflned aquifer 
la thought of to be concentrated at the bottom of the aquifer. 

The unit of resistance to flow is same as that of time e.g 
days 

Alike the inter-aquifer flow, the flow between the unconfined 
aquifer and the neighbouring stream experiences resistance if 
the bank of stream forming boundary is clogged This resistance 
is given by 
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(3 22 ) 


C = 

A 


/ k. 


C « 

B 


b 

2 


/ k 


(3 23) 


3.B PIEZOMETRIC HEAD IH UNDERLYING AQUIFER 

The plezometric head at any point In the underlying aquifer 

ia the height to which the water will riae if a piezometer ia 

introduced at that point in the aquifer. The piezometric litie 

conaiata of a number of pointa repreaenting the piezometric 

head in the underlying aquifer. In fact the aquifer will be 

connected at ita enda to teaervoira or water bodiea located 

uaually at far diatancea on either aide Under the ateady state 

condition the piezometric head line will be given by 


h -h 

h = — I — X + h, 
^ L 

I 


(3,24) 


where 

h^= water surface elevation in the left reservlor.L; 
h = water surface elevation in the right reservior,Lj 
L = distance between the two res ervoira , L . 

In ease the flow ia unsteady due to the fluctuation of water 
levels in the bounding reaervoire h^(x,t) can atill be found. 
But under practical clrcuatances such data are seldom known. 
Presence of tubewella etc. create additional complexities. Since 
we are intrested in a small region, the piezometric surface for 
practical purposes is well given by 


h = -±iU= X + h , . . (3.25) 

o 1 «*». 

X 

where 

h piezometric head at right end,L; 
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h . plezomfetrlc head at left end , L 
This can be written as 

h = px+q (3 26 ) 

O 

Describing h as function of x only is not a limitation of 
the method presented in this study for finding a aolultion for 
the phreatic surface Transient piezometric surface can also be 
treated The only limitation la that it should be a known func- 
tion of space and time. This means that it should not be affec- 
ted appreciably due to the int et -aqui £ er flow In practice the 
presence of a seperating layer of low permeability makes such 
ati assumption valid 

3.6 TRUST ENT UNSATURATED FLOW ABOVE THE WATER TABLE 

The recharge to water table from top will take place either 
due to rainfall or excess irrigation or artificial recharge of 
water. This recharge may be assumed uniform but will be unsteady 
usually The recharge is assumed steady in steady state problems 
or for the sake of simplification in case of unsteady state 
problems, but it should be noted that even a steady rainfall 
will cause an unsteady recharge to the water table due to the 
action of unsaturated zone The presence of plant roots will 
further complicate the problem. Negative recharge l,e evapora- 
tion from water table will take place during interstorm periods 
which will again depend upon a number of meterologlcal , hydro- 
logical and topographical parameters. 

3.6.1 GOVERNING DIFFERENTIAL EQUATION FOR UNSATURATED FLOW 

The complex process can be fairly approximated by a partial 
differential equation, usually called Richard’s equation with an 
additional sink terra added to it to represent the transpirat ion 
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loss The equation can be derived easily by combinine Darcy’s 


law with continuity aquation as follows 
According to darcy'e law 

a = - k ^ 

The unaaturated head la given by 


(3 2 7 ) 


(3.28) 


where 

q^ rate of unsaturated flow in positive z-dlreetlon LT ^ 
h piezometric head in unsaturated soil.L, 

U 

unsaturated permeabl 1 i ty , LX” * , 

z vertical co-ordinate measured from surface vertically 
downwarda.L , 

p^ soil matrix potential ,negativeve in unsaturated soils, L 


From continuity, 


where 




az 


+ T + 
a 






0 


(3 29) 


& moisture content of soil mass , dlmens ionl ess , 

-1 

T rate of transpiration per unit depth, T; 

a 

t time,! 

Combining the two, we get 




or, fk — ^ — [ “ “ T “ — • (3.31) 

The governing partial differential equation becomes 


a ^ ah ^ ae 

- — k — + — = 0 

ae 


(3 30) 


ak. 


a& 


26 



s 


SQ 


(3.3e> 


6 - 1 - _H - T - - 

&z ^ Sz ® C>t 


Sz 


whete i5 IS called diffusivity of the soil winch is the rate of 


flow du« Lo arid gLv«n by 


dp 


U 


S = \ 


az 


C3.33) 


Thus the first term on left hand side represents flow due to 
soil matfix potential where as the second terrri represents flow 
due to gravity and the third term the loss of water per unit 
depth due to transpiration by plants. The termi on righ hand side 
represents rate of change in storage. The rate of recharge to 
water table will be given by 


w = - <5 ^ 


z=l 


( 3 . 3 ^ ) 


where 1^ is the depth of water table from surface. The recharge 
to water table can be evaluated if the above equation can be 

solved. 

3.6.2 ASSUMPTIONS INVOLVED IN DEPRIVING THE UNSATURATED FLOW 
EQUATION 

The above derrivation makes the following assucriptions t 

(1) Darchy's law is valid for unsaturated flow. 

(2) Isothermal conditions exists in unsaturated flow sone. 

(3) t- ^and p^are single valued functions of moisture content d. 

(4) Incompressible flow takes place through a nondef ormable 
media . 

(3) The air effects are negligible. 

(6) The roots are distributed uniformly in the root zone. 
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Fie 3 3 Definition sketch for unsaturated 
ground water flow 






363 INITIAL AND BOUNDARY CONDITIONS FOR UNSATURATED FLOU 
The boundary conditions at the top and at the water table ae 
well as the initial condition at e as iollows 
3631 Initial Condition 

The initial condition in its most general form may be written 
as 

a = f(z) at t = 0 - (3 35) 

3 6 3.2 Boundary Condition At The Water Table 

The boundary condition at the water table is always given by 

& = & at'TL =1 ta 0 

e z, 

(3 36) 

3.6,3 3 Boundary Condition At The Ground Surface 

The boundary condition at the ground surface will vary with 
time and to define it, the entire process is to be understood 
When the rainfall starts, initially the entire rainfall gets 
infilterated till the dry soil surface gets saturated. Then the 
rainfall in excess of infiltration gets accumulated on surface 
filling up the depressions and then run-off starts The run-off 
continues to take place till a short time after rainfall stops 
Then evapouratlon from the soil surface will consume the volume 
in detention storage alongwith infiltration During all this the 
soil surface remains saturated Thenafter desaturating of soil 
will occur due to the evapouratlon till the next rainfall takes 
place . 

The boundary condition during the first phase prior to soil 
becomes saturated is given by 
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j; X I 

-6 - 5 — + k 

&z u 


tv ^ 


t < t 


Bi 


(3 37 ) 


During ih* second phase soil surface remains saturated and 
the boundary condition is given by 

d = e t:^t<t. C3 38) 

a SI di ' 

In the third phase the rate of evaporation from soil surface 
can be found from met ereological data and then the boundary 
condition is gLv«n by 


6 


d-e 

Sz 



t . < t 

di 


'be 


(3 39) 


where 


r transient rainfall rate.LT *; 

e potential evaporation rate,LT”* , 

t, time Instant correspond Ing to begining of first rainfall, 

Di 

t^2 time instant cot responding to begining of next rainfall, 

t time instant corresponding to end of first rainfall; 

t time instant corresponding to ®r,d of next rainfall; 

62 

t time Instant corr espondlng to the time when surface gets 

saturated during first rainfall, 

tdi time instant corresponding to the time when the surface 

starts drying after first rainfall, 

t duration of first rainfall; 

ri 

t duration of one rainfall period which is the time since 

Ci 

the begining of rainfall to the beginning of next rain; 
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In case the rainfall la not auffeclent to saturate the aoll 

surface, the phase I will change directly to III and phase II 

will be absent during that rainfall period The phase I boundary 

condition will prevail till t^^ and the phase III bounda*'y 

condition will continue till t „ Value of t will be obtained 

e2 Bi. 

during the solution but t^^ requires to be defined as tha tnput 
data If volume in detention storage is assumed to be zero. 


Otherwise t^ can be found from following 
assumed that the volume In detention storage 
potential rate, which is given by 


(3 40 ) 

equation, if it is 
evaporates at the 


V 


d 


6 


t 1 

di r4 



(3 41 ) 


where i la the rate of infiltration given by 

^ " H n * “u ] z=<. 

This requires knowledge of detention volume ( expressed as 
depth of water) and rate of evaporation under prevailing condi- 
tion Later one is also required for defining boundary condition 
during third phase. 

364 SOIL MOISTURE CHARACTERISTICS 

The unsaturated soil moisture pressure p^, unaaturated permeabi- 
lity k as well as dlffusivity 6 are all nonlinear functions of 
u 

moisture content & and hysterisis la exhibited. In the present 
study hysterisis will be neglected and it is proposed to use 
the average S ~ 6 , , and k^- B curves. 

Various investigators have investigated the p^ - © and k^ - 6 

relationships and proposed a power form relationships as shown 
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in the table 3 1 


TABLE 3.1 

SOIL MOISTURE CHARACTERISTICS 


REFERENCES 

RELATIONSHIPS 


3 5 

Averjanov (1950) 

11 

Iramy (1959) 

K = 


-C2+3X) 

Brooks and corey (1966) 

„ r »u 1 


■t-J 


-X 

^ = 


o ois w -i- a o 

Mualem (1976) K = ( ^ ) 

W = J V, ds 



Mualem’s (1976) formula i.a baaed on generalised macroscopic 
approach and la found to be applicable for a wide variety of 
eoila and tnolature contenta greater than wilting point It 
employes a macroecopic physical property which indicates the 
amount of work per unit volume of soil mass required to drain a 
saturated soil to the wilting point 
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Eagleson (1976 a) indicates that the power function model 
by Brooks and Corey (1966) applies nicely to the data of Talsma 
(1979) and Moore (1939) Also he indicates that a relationship 
between index of in K - relationship and \ , the pore size 
distribution index la valid for a wide variety of soils From 
the discussion it can be taken as 

K = K^(<^)^2 (3 43) 

P = P^(?!.)^2 (3 44) 

Therefore the diffualvity of the soil can also be represented 
by a power function relationship. 

A = Dj^(^!-)^2 (3.45) 

where 

K non-dimensional unsaturated permeabi li t y , 

P non-dimensional pressure head,p^ / p^ 

A non-dimensional diffualvity, 

bubbling pressure 

365 ACTUAL TRANSPIRATION RATE 

Corroborated results by Denmead and Shaw (1962) suggest that 
for a given potential rate there is a threshold average coot 
sone soil moisture content ( or the corresponding average soil 
moisture auction ) below which the transpiration rate through 
the plant falls below the potential as indicated by Cordova and 
Braes(19Bl) . 

Neghasi (1974) suggested following relationship to relate the 
transpiration rate with soil moisture content, 

t^ / tp = ^ (3 46) 

t^ / tp = 1 4^ ^ . (3,47) 
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where 


Actual tranepi rat ion rate 
Potential transpiration rate 


t 

a 


* 

^ Threshold moisture content 


Eagleson (197Sb) describes a transpiration model as 


^a / 

t = (P-P*)/r(r -V )y T 1 

p ' L 8 (ji pvj 

P < P 

(3 

48) 

"a / 

t = K 

P V 

A 

P * P 

(3 

49) 


where 

r resistance to moisture flow in soils, 

r resistance to moisture flow in plants, 

P 

Y, specific welfiht of soil; 

W 

* 

P critical leaf moisture potential, 

Ky transpiration efficiency 


The values of various coefficients can be found from plant 
physiological literature including the works of Evans (1963), 
Cowan (1968) Slatyer (1 9 6 9 ) and Kramer (1969) 

The transpiration model used in this study is given by 

t^ / tp = e'j^ + e'^ 

t / t = 1 a (3 51) 

a p 

3.7 MUMERICAL METHOD USED FOR VARI FI CATION OF RESULTS 


Predictor - corrector method of Douglas and Jones (1963) has 
been used in this study to find numerical solutions, for both 
saturated as well as unsaturated flow . 

The partial differential equation governing the saturated 
flow in nondiment ional form is given by 

d H a n 2 a H 


^ T 


H 5 — + ( ) + I 

d X J X d X 


C H + C H + U 
o o o 
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. .(3.52) 



where H,X,T,S,C^ and U are non-dimensional quantities represen- 
ting height o£ pheratlc sur £ ace , lior izont 1 distance , time , slope, 
vertical resistance and recharge respectively. The quantities 
with respect to which nondimensional i sat ion has been carried 

out depend on type o£ problem and hence will be discussed when 
a particular problem will be considered 

Similarily, the equation governing unsaturated vertical £low 
in nondimensional form can be written as 

& 4, 

= A a'( - k' . (3 53) 


Both the equations can be expressed as 




d X 


= ?( X,T,U,U' 3 


^ T 


(3 54) 


where U is the dependent variable, X Is the independent variable 
and T is nondimensional time In both the cases ’F’ can be 
written as 


F = gj^(X,T,U,U* )■■ ■ + g^CX.T.U.U') (3 55) 

which is same as eqn.(1.9) of Douglas and Jones(1963). 

Using Douglas Jones Predictor-Corrector method , the predictor 
is given by 






: T u 

n+*yii' 




1 


0.5(U^ 





(3.56) 
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The corrector la given by 



Thua a eet of linear alntul taneoua equationa are required to 
be aolved for predicting the valuea of U, ^ from vhich the 
valuea at the end of time atep can be found by aolving another 
set of llnar aimultaneoua equations The value of H In case of 
steady state problems has been obtained as T oo 
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CHAPTER 4 


SOLUTION DISCUSSION AMD APPLICATIONS 


4. 1 EXACT SOLUTION FOR STEADY STATE FLOW OVER HORIZONTAL 

SEMI -PERVIOUS BED 

Leakage to or from the top moat aquifer (we are Inteteated 
in the flow in the top unconfined aquifer) may take place in 
addtion to recharge from top whenever the underlying bed ie not 
impervious The quantity of artisan flow will be proportional to 
the difference of heads in utiderlylng and unconfined aquifer 
Since the pheratic surface is a curved one,artiaen leakage will 
not be uniform and hence cannot be treated as recharge from top 
which is usually taken as uniform Noticing this fact Uesseling 
and Uesseling (1986) have found a linearised solution An exact 
solution can be found using the solutions provided by Sikkema 
and Van Dam (1982), and Nieuwaal and Nijhuis (1974). 


4.1.1 NON-DIMENSIONALISATION 

The governing partial differential equation for steady state 
flow over horizontal semi-pervious bed Is 


|_(k h|^) + — ) + ie = 0 

dx s dx c 


C4.1.1) 


Using the non-dimensional parameters, 


H=h/|h 


(4 1 2) 


X=x/L 


(4 1.8) 



2 


c = ** 

o K CTTT 

a A ’ 


£=1/C 

' O 

where is the equivalent artesian head given by 

h = (h +UC ) 
e ' o ' 


(4 1 4) 
(4 1.5) 

(4 1 6) 


eqn (4 1 becomes 
JL.C H £5 

«»x I ” dx J 


C H + C 
o - o 


(4 1 7) 


(Use + when h >0 and - when h <0) 
' e e 


The solution o£ the above equation is found for following 
cases 

(i) Uhen heads at both boundaries are specified , 

(ii) Uhen head and flux at one boundary are specified 
412 SOLUTION 

The eqn (4 1.7) can be written as 


HH ^ (HH ) = C H ( H + 1 ) 


(4 1 ft) 


dH 


(HH* )^ = I C^(H^ + 1 5 H^+ a) . (4 1.9) 


HH 




= A(H) 


. (4.1 . 10) 


(Use -ve sign for H >0 and +ve sign for H <0) 


where 




(4.1 11) 


and g(H) is given by 


g(H) = + 1.5H^+ a 


(4.1.12) 
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The solution o£ eqn (4 1 9)depends on the value of 'a' which 
is governed by the relative elevations of artesian head line 
and water levels in channels forniing boundary 
4 12 1 Uheii both the boundaries are specified 

CASE-I WHEN = 1 

In this ease H = 1 will be a solution of eqn (4 1 S) which 

means water table will be horizontal one with its height equal 

to h +w c .Fig (4 1 1) shows the profile, 
o 

Before considering other three eases for h >0, h =0 and h <0 
it is worthy considering flow in semi-infinite aquifers 
1 2 a,b,c,d) 

When h >0 (/f<f-4.1.2 a,b)the height of phreatic surface changes 

A 

from h^to continuously and it meets the h^-llne tangentially 
at certain distance from boutidary Boundary effects are confined 
In a distance 1 only and beyond that water table is horizontal 

oft 

Selecting origin such that 


(4.1 13 a) 
. (4.1.13 b) 

From eqn (4 1 9) with eqn. (4.1 13 a) we get 

a=+0 5 (4 1 14) 

The solution for this case follows from Nleuwaal and Nljhuls 
(1974) as 


X=0 


H=1 




ax 


= 0 


X=L , H=H 

eS B 
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Flfi. 4.1.2 (b) 


Fig. 4.1.2 Flow ovet^ afeini-pebvious bed in aemi-inf inlte domain 
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Fife. 4 1.2 (d) 


Fig. 4.1 2(d) Flow over aeai-pervioua bed in aeal-inf inlte domain 
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VH+O 5 


+ •/1/6 In 


•/H + 0~5 — 


y"H + 0“ 5 + V i 5 


/X7 


X + b 


(4 1 15) 


where b ia a constant of integration Since h— >h asymptotically 
b is taken from x=0 at h=0 99h or 1.01 h Value of 'b' thus 

« A 

comes out to be 1 39 Now the profile can be computed from 
equ (4 1 15) and 1 is given by 

CS 


1 = / 6k c h 

C* yf * * 


VK +0.5 + 


In 


■/H +0 5 - -/l~S 


yir+o 5 + vTTT 


+ 1 39 


. C4 1 16) 

Uhen h <0 (/ttf-4 1 2 d) the height of phreatic surface changes 
e 

continuously from h to **ro as it meets the horizontal bed at 
certain distance from the bousidary Selecting origin such that 

X=0 , H=0 <4.1 17 a) 


X=L 


oB 


H=H 


(4 1.17 b) 


From 6qn.(4 1 9) and eqn.(4 1.17 a) we get 

a=0 . (4 1 IS) 

The solution will be a parabola given by 

r 1 ^ 

H = g- X + V 9/C^ - 1 5 ..(4.1.19) 

where the second constant of integration 'b' has been evaluated 
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frojn^eqn (4 1 17 a) The profile can be found frofli eqn (4 1 19) 
and L ^2 given by 

ff 

L = ■/6k ch yi 5 + H - Vl 5 (4 1 20) 

s e I ^ 

J V (t. 

Uhen he = 0 , non-dlmensi onal laatl on with respect to will 
not be possible In fact in this case, using the dimentional 
variable itself and integrating equation 

W .2 

h h + h - h/c = 0 (4 1 21) 

we get, 


i»h 

( h — ) = 
&x 



(h^+ a) 


(4 1 22) 


In this case (/itf-4 12c) 
changes from h at boundary to 

ft 

bed tangentially Selecting the 

x=0 , h=0 

. h=h2 

From eqn (4 1 9) and eqn (4 1 

a=0 


the height of phreatic surface 

zero and it meets the horizontal 
origin such that, 

, dh/i»x=0 . (4 1 23 a) 

. . (4 1 23 b) 

:3 a) we get, 

C4 1 24) 


The profile in this case will be a parabola given by 


h = 



(4.1 25) 


The second integration constant comes out to be zero Profile is 
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by eqti ^4 1 24) fexplldltly and L is given by 

CB 


1 

C2 



(4 1 26) 


CASE- 1 1 WHEN h > 0 
e 

L and L are first found from eqn (4 1.16) for h >0 
Ci ca e 

If 1 i the profile can be computed from eqn (4 1 15) 

Otherwise following cases arise 
Case- 1 1 A When l>Hj^>H 2 

In this case (yftf-4,1 3a) maximum water table height h < h 

me 

will occur somewhere in the flow domain. Selecting origin such 


that 


X=0 


H=H 


dv. 


m 


= 0 


X=L 


H=H 


X=L . H=H 


* ' t ' a ' * 

From eqn (419) and eqn (4 1 27 a) we get 

2 


from which 


a = -H + 1 5H 
m m 

0<a<0 5 


and therefore equation gCH) = 0 can be written as 

g(H) = CH-rJ(H-r^ )(H-r^ ) 

Since 0<Hm<l , 0<H<Hm , the solution is given by 



X + b 


where 


8ln> = (H-r^)/(H^-rj) 

2 

sin a = (r 2 -rj^)/ (r^-r j) 
on ^ 


(4.1 27 a) 
(4 1 27 b) 

(4 1 26) 


(4.1.29) 


(4 1 30) 


(4 1.31) 
(4 1 32) 
(4 1 33) 
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Fig. 4.3.1 (b) 

Fig 4 13 Flow over horizontal etemi-pervious bed due to a 
positive artesian head 
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Cc) 


Fig C4 1.3) Flow over a horizontal eemi-pervious 
bed due to a positive artesian head 


Fig 413 


(contd.)Flow over horizontal 
positive artesian head 


sejni-pervious bed due to a 
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and ’b’ Is an arbitrary constant of integration which can be 
evaluated by using equation (4 1 27 a) so that 

b = Cr^-rj)E(90° ,«) + r^F(90'^,£A) (4.1 34) 


E and F 


are elliptical integrals given by 



d& 


dd 


(4 1 35) 


(4 1 36) 


Case II B Uhen Hj>l>H2>0 

In this case (.fig-k 1 3b) the pheratic line will meet h line 

© 

somewhere but not tangentially Selecting origin such that 


X = 0 , H=1 , = 0 

i»X 

X=L^ H=H^ , X=L2^ . H=H^ 

From eqn.(4 1 9) and eqn.(4 1 37 a) ve get, 


ah 




-05 


. (4 1.37 a) 

(4 1.37 b) 


(4.1 38) 


ftom whlcli a > 0 5 

and therefore equation g(h)=0 can be written as 

gCh) = (h-r^) (h^ + 2ph + q) . . (4.1 39) 

2 

Since 00 > > 1 and 1 > > 0 and noting (q > P ) 

the solution is given by following equation for either part 
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- 2 yr 


+ 2 

where 


and b is 
eqn (4 1 

b = - 2 


where 


E C«^ M ) + 



] FC<^ / ot) + 


Sin / 1- 

= ± 

(1 + Cob 


A + r - h 

Cos ip = * (4. 

A - r + h 





r > 0 

i 

I constant of integration which can be evaluated 
37 a) and eqn (4 1 40) as 

1 

■/IT E(4,^ ,«) + (VTT + ) FCip .a) + 

* v~ir * 


« « A - r - p 

Sin^ot = = ‘ 

2 A 


C4 


C4 


+ 


Sin 
2 tTA — 


ip^ / 1- K^Sin^^ 
( 1 + Cos ) 


is given by 


COfi 



A + r - 1 

A 


A - r + 1 

i 


(4 


(4. 


1 40) 


1 41) 


42) 


1 43) 


from 


1 44) 


1.45) 
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Cfiif.e II C Uhen H >K >1 

4 g 

In this case 13 c) minimum water table height h^>h^ 

will occur somewhere in the flow domain Selecting origin such 
that 


X=0 


H=H 


^h 

5x 


= 0 


X=*L. 


H=H 


X = L 


2 * 1 ’ 2 
From eqn.(4 1 9) and eqn (4 1 46 a), 


H=H 


(4 1.46 a) 

(4.1 46 b) 


2 


a = H ( 1.5 - H ) 
m m'^ 

Uhen H =1 5 In this case a « 0 


and therefore the solution will be a parabola given by 

C 


H 



5 


(4 1 47) 


(4 1 48) 


Uhen H < 1 5 In this case 0 5 > a > 0 

m 


so that gCH)=0 can be written as 

g(H) = (H-r^)(H-r^)(H-r^) . (4.1 49) 

and therefore the solution is given by 


r^ F(^ /oi) + (r^- tj^) 


E(^ / a) + tan 


where 






X + b 


(4.1 50) 
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Sin « 


1 


(4 1.51) 


r - r 
a 2 


r ~ r 
a i 


r = H 

a m 


(4 1 52) 


tan ip 


/!ZZ 


*^3 " *^2 


(4 1 53) 


and the second constant of Integration 'b' can be evaluated 
front eqn (4.1 46 a) as 

b = 0 ..(4154) 


Uhen H >15 In this case 
nt 


a < 0 


so that equation g(H) = 0 iho-b fine--' ..-ryodt 


and therefore the t .lution la given by 


-!■/ 


A E(<^ / a) + 


•/k +- 


■/IT 




) + 2yA ■/ 1-K^Sin^ 




1 + Cos ip 


/ 


2 C 


X + b 


(4.1 55) 


where 


A a 


c + 2<i 

i. 


. (4.1.56) 


CV 


/U 


» . 


A + - h 


A “ r^ + h 


(4 1.57) 
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K 


2 


Sin ^ 


A - r - p 

L 

~5a 


(4 1 58) 


r 


1 



(4.1 59) 


CASE-III UHEN h < 0 
= e 

L and L are first found from eqn (4.1 20) for h < 0 

C 4 02 e 

If 1 i the profile can be computed from eqn (4 1 19) 

Otherwise fflinlmum watei table height h <h will occur somewhere 

me 

in the flow domain Selecting origin such that 

X=0 , H=H , dH/dX=0 (4.1 60 a) 

m 

X = L^ , H =H, , X = L, , H=il, j 

From eqn (4.1 9) and eqn (4 1 60 a), 


2 

a = H ( 1 5 - H ) 
m m 


(4 1 61) 


Uhen H =0 5 In this case (/itf-4 14a) a = -0 5 

m 

so that equation g(H)=0 can be written as 

g(H)=(H+l) (H-0.5) (4.1 62) 

and therefore the solution is given by 


y H-0 5 - y “2 /T ■ tan"* -/Z/S H-l/l" = Vc^ /6 X + b. (4.1 63) 

where 'b' la a constant of integration which can be evaluated 
by using eqn (4 1 60 a) as 

b = 0 . (4.1.64) 

Uhen H >0 5 In this case (/f^-4.1.4 b) a < -0 . 5 

m 

so that equation g(H)=0 can be written as 
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Fie. 4. 


t 



Fifl,4 1,4 (a) 



Fig. 4.1.4 (b) 


.4 Flow over horizontal aefni-pervioua bed due to a 
negative artesian head 
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Fie. 4.1.4 (c) 


Fie 4.1.4 (eontd ) Flow over horizontal a emi-pervi oue bed due to a 
negative arteaian head 
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C4 1 6b) 


e(H) = CH-r^) (/+pH+q) 


and therefore the solution ia given as 



E /a ) + 



F(4> / a) + 


Sin ^ / 1- K^Sln^^ 

+ 2 -rr- = ± 

(1 + Cos iP) 

wher 6 

A + r - W 

Cos <p = ^ 

A - r + H 


X . b 

y 3 


A 



+ 2 q 


Sin « = K = * 

2 A 

r > 0 


and 'b' is a constant of integration which can be 
from eqn.(4.1 66) and eqn (4 1 60 a) as 

b = 0 

UHEN H < 0.5 In this case (/f^-4.1.4) -0 5< a < 

in 

ao that gCH) can be written as 

g(H) = (K-rj^)(H-r2)(H-r3) 

and “1.5<r^<-l , -l<r2<0 and 0<r3<0.5 . Since H>l,the 
is given by 


(4 1 66) 

(4 1 67) 

.(4.1 68) 

. (4.1.69) 

(4.1.70) 
evaluat ed 

(4 1 71) 

0 

(4 1 72) 
solution 
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2 2 
1-K Bin 4> 



X + b 


(4 1 73) 


where 


K 


2 


1 




= ein 


2 


oi 


tan ^ 



(4 1 74) 


C4 1 75) 


r. = (41 76) 

3 m 


The second constant of Integration 'b' comes out to be zero 

CASE-IV UHEN h = 0 
e 

1 and 1 are first found from eqn (4 1 26) for h >0 

CA Cft ® 

I f 1 5: 1 +1 the profile can be computed from eqn (4.1.25). 

X c* ca 

Otherwise following cases arise • 


CASE-IV A UHEN h >0 , h >0 

4 2 

In this case (/itf“4 15a) minimum water table height h^ > 0 
will occur somewhere in the flow domain. Selecting origin such 
that 

_ . . dh 

■ ■> ' -S5r “ “ (« 1 77 .) 

x=l^ , h=h^ , *"^8 ' 1-77 b) 

From eqn (4 1.22) and eqn-(4.1.77 a), 

a = -h^ C4 1 78) 

m 

so that equation g(H)=0 can be written as 





fi(h) = (h-hn») (h^+h£nh+h^ ) .. (4 1 79) 

in 

and hence the solution is given by, 

-2 ( 3 / E + 1^3)" + ' CJJ ^ ^ 

\ ^ 

( 

s > 

* 2 Wl = 1 " 

- - — -- - - . ■_ ■ 

(1 + coa <f>) 



X + b .... (4 1 00) 


where 


cx = 15 


. . . (4 1.81) 


h +h -h 
V y mm 

coe<p = j' *" 

/JS h -h -h 
y mm 


. . . (4.1.82) 

and 'b’ is a constant of integration which comes out to be 


b = 0 ... (4 1.83) 

CASE- IV B 

In this case dh/dx WJill not be zero anywhere in the flow 
domain since 1 < 1 .Selecting origin such that 

x=0 , h“0 ..(4.1.84 a) 

x=l , h=h ..(4.1.84 b) 

X 2 

From eqn (4.1 22) and eqn.(4.1 84 a) a=0 

and therefore the solution will be 

2 

h =-^ (4.1.85) 

From above it is seen that dh/dx^O at x=0 which is not 

possible . This is because such a situation cannot occur as there 

will be seepage face at the other end 

B7 





Moat of the aolutlon found above Involve independent variable x 

or X expressed as function of dependent variable h or h/|h | 

0 

Uhen heada at both the boundaries are apecified, a trial and 

error procedure la iieceaaary to find h Once h la kirowii , tlie 

mm’ 

profile can be easily computed 

4122 Uhen Head and Flux are specified at one Boundary 

Uhen Head and Flux at one boundary are known, procedure Is 

d h 

straight forward First constant a is found using h and , 

(hh* -h^ +1 5h .(41 86) 

e 

Now depending on h > 0, h < 0, h = 0, and depending on the 

6 0S 

the values of h^ ,h^ and a the proper case to which the problem 
belongs can be found from the last section From that solution, 
the profile can be plotted from the reference point to the 
boundary at which h = h^ , let this distance be 1 Then the 
distance to the other boundary will be (1^" 1^) where 1^ is the 
length of entire domain Putting x = 1 - l^ln the solution will 
give us the value of h. , thus profile at other side can also be 
plotted 

4.3 1 APPLICATION TO DRAINAGE PROBLEM 

4.1.3 1 Design and Analysis problem 

Usually drains are employed to take away the artesian water 

whenever artesian head is high . This corresponds to case II A 

as H >0. 
e 

In a design problem, spacing and depth of drains are to be 
found to maintain maximum water table at certain level . From 
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equation 


'x' / ‘ •'e ■'e ■= 'n (4 1.87) 

where I la given by 
n 

/ o() - bj /6 /Cr^-r^) 

. (4.1.8S) 

r, ij!» and c* follows from eqn (4 1 29) to eqn (4.1 31) Assuming 
a suitable depth of drains 'd* , spacing of drains 1^ can be 

found for a given h Thua 1 can be calculated for various 

m X 

values of d and a proper combination la selected The criterion 

for this selection may be economic In an analysis problem the 

maximum height of water table and water table profile is to be 

found, knowing ’d’ the depth of drains and spacing • 

The procedure will be a trial and error one and is as follows : 

1. Assume h and find h /h and d/h . 

m m' e ' e 

2 Find spacing and check whether it is correct or not If it 
is not correct, try another value of h Once h is found, rest 
of the profile can be easily computed 
4132 Comparison to Linearised Solution 

Uesaeling and Uesaeling (1986) gave a linearised solution of 
eqn. (4 1.1) which in our notation after some rearrangement can 
be written as 

H =1+(D-1)/ cosh (05/ Jc D ) (4 1.88) 

iT\ ' * y 

FIG (4 1.6) shows H calculated from linearised and exact 

m 

solution for different depth of drains d/ The dotted line 

shows curves for linearised solution and dashed lines for exact 
solution. The errors are plotted against (/itf-4.1 7). 

Errors are large for small seen that linearisation 

results in overestimaLion of h^ and this overeat imat ion is 
larger for smaller values of * It is found that error in 


60 



found from linearised solution is less than 5% for D 0.5 and 

£^^02 Also for D ^ 0 2 maximum error in h is only 5.47% 

m 

for £ > 4 . 

Thus it can be concluded that linearised solution works well. 
Since linearised solution is simpler one, it should be adopted 
for use Exact solution should be used when c is less than 4 and 
only when 0^05. 


4 1 3.3 Comparison with Young's Solution 

Due to Dupuit's assumpt ions , the solution based upon it fails 
to account for convergence near a partially penetrating drain 
and seepage face youngs(19fid) found a solution which takes into 
account these factors The solution can be well described by 
following empirical formula 


-1 65 L / 2 h 

= D + (1 - D)(l - 6 ° (4 1 89) 

Youngs' formula does not take into account rainfall recliarge 
ft om top and resistance of underlying stratum which in practice 
is found to occur usually Solution for this case using Dupuit’s 
approach will be 

To account for convergence , Hooghoudt’a equivalent depth 
'D ’ should be used Instead of 'D' This is computed by using 


eq 

simple formula for Hooghoudt’s equivalent depth given 
and Antonopoulos (1981) which are as follows . 

D = — ^ 

eq 8^ 


by Sakkas 


(4 1 91) 
(4 1 92) 
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-1 


(4 1 93) 


o 

<4*^= 0 733 loe<y - 0.354 (4.1.94 a) 

0 01<<y<2 


= 0 . 345 ( loeo‘) + 0 . 508 logo* -0.326 

(4 1 94 b) 


vhaf 6 


Cf 


0 1250- + 0 .160 
ct 


L 

D h 


e 


-0 443 


(4 1 94 c) 
(4.1.95) 


Using Young's formula, and exact and linearised formula without 

and in combination with Hooghoudt’s equivalent depth, L/h^ and 

h„/h for various values of D/h were calculated It was 

found that linearised solution in combination with Hooghoudt’s 

equivalent depth compares with Young’s solution nicely Assuming 

Youngs' solution as correct one, the errors in the linearised 

solution In combination with Hooghoudt’s equivalent depth ate 

in getieral leas than 5 % though maximum error la 6.5 % in a 

particular case PIG. (4 1 8) shows the plot between h^/h^ vs. 

L/h for various d/h for these two solutions Exact formula 
o o 

with Hooghoudt's correction gives better result when o = 0 2 as 

maximum error la only 3 26% at L/h^ =05. It is noted that at 

d/h 0.1 the Dupult - Forchmeir approach fails as flow will 
o 

become more and more vertical Therefore it was concluded that 
linearised solution Jn combination with Hooghoudt’s correction 
should be adopted for use This is because this formula can be 
used for a wider varlely of problems when rainfall la there as 
well as when a aemipermeabl e stratum separates the upper and 
lower aquifer which is commonly found to occur Jn practice. 
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0 . 12 

0.09 

0.20 
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0.3 
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mm 
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B 
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B 
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B 
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2 . 30 

0 . 98 

B 

B 

- 0.19 

- 0.2 

- 0.20 

0 . 7 

B 

B 

B 

2 . 36 

1. 24 

0-50 

0 . IG 

0 . 13 

0 - 09 

0 . 12 

0 . B 



2.01 

B 


0.47 

0 - 17 
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0 . 18 

0.9 

- i , 16 

0.91 

0.95 
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4. a LINEARISED S»LUTIOKI FOR STTEADY STATE FLOW OVER SLOPING 


SEMI -PERVIOUS BED 


The steady state flow over a eloping e eml -pervloua layer 
haa been considered by Nualem and Bear(197S) .They derived 
an anlytlcal aolution by transforming the partial differential 
equation and then linaerlalng it. Their solution was only for a 
special case when artesian head was equal to ,the thickness 

of sloping aemlpervloua stratum In the present case arbitrary 
artesian heads are considered. The artesian head line will be 
assumed parallel to bed for flow between parallel drains and a 
linearised aolution will be found The aolution will be compared 
with the numerical solution to examine its validity. 


4 2.1 NON-DIMENSIONALISATION 

The partial differential equation applicable in this case 
will be 


J (k h 
s 

Sx 





I dh 


dx 





u 

+■ u> 3-0 


h(0) = d 
h(L) =d + I 

Defining effectlce artesian head as 

h = h + c w 
e o 

and using non-dimensional parameters 


H = h / h 
' e 


X 


S 


X / 1 
K c h 

3 < 


/I 


2 

X 


I L / h^ 


. C4.2.1) 

(4 2.1 a) 
. (4.2.1 a) 


(4.2.2) 


(4,2 3) 
(4 2 4) 
(4 2.b) 
(4 2.6) 
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V6 get 


a 


c H + S dH 

ix“ 


-C H + C s 0 
o o 


(4 2.7) 


with, 


H(0)“D (4 2 7 a) 

H(1)=D+S (4.2 7 b) 

The exact aolutlon of eqn (4.2.6) la difficult to find. 


Hence 1 ineariaat ion approach becotnea inevitable 
4 2 LINEARISATION AND TRANSFORMATION 

Eqn (4 2 6) can be lineariaed to get, 

D H' ’ + H’ -CqH + 0 

Using the trana format ion , 

H-U e’”* 

we get , 

DU" /n - Cq e"*^ /n =0 
U(0) = D 


. . (4 2 8 ) 


. (4.2.9) 


. (4 2 10) 

(4 2 10 a) 


wher e 
and 


U(l) 


(S+D) 


(4 2 10 b) 


m = S/2D 
n = m D ' Cp 


(4 2.11) 
(4 2 12) 


423 SOLUTION 

The solution of the above eqn (4 2,18) can be written aa 
U = H Sinh ^ n X~ + Co-gh 

(4 2 13) 


S = D-1 


(4 2.14) 
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CD-I) Cosh 


(4 2 15) 


atid 




m 


CD + S-1) e - 


y/ n/D 


Slnh Sy n/D 


The maximum height of water table will occur where 

& H 


0 


C4 2 16) 


S X 


from which the distance to maximum height of water table X 

can be found to be given by the root of 

m C - C / 

a 4 -/ D 


m 


tanh 


X J 
my D 


= 0 


C4 2 17) 


m C 


- c r^ 

4 ^2/ D 


and then H will be given by 
m 

*<a ^TKTT * Cj Cosh , ^-t.X ^ ^ 

. c^ 2 la) 

4.2 3 COMPARISON WITH NUMERICAL SOLUTION 

The maximum height of the water table h^ and its location 

are computed from the eqn.C4.2 26) and eqn.C^ 2 2 7 ) r eapect ively 

and also by a numerical approach for a wide range of parameter 

a varying from 0 2-10 The computation was carried out for five 

values of non-dlment ional slope S The error In h^the linearised 

solution was found and presented in tables From there it is 

seen that errors in H are small as long as D ^0.4 Errors in 

m 

X are small but increase with n and when ground water level 
m 

reaches near upper drain errors become unacceptably large. Of 
course such situation will not be found in field as then the 
spacing will be uneconomicaly small. Thus for field situations 
linearised solutions are found to be acceptible 
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TAbLE 4.2.1 

ERRORS DUE TO LINEARISATION FOR S =0.05 



OJ 

m 

O 

' 0.4 

o 

■ 

0.8 

1 

2 

4 

6 

a 

10 

0.1 

43 



68.2 

63.2 

36.2 

5.39 

- 27.3 

- 16.9 

- 21 . 7 ^ 

0.2 

26.9 

31 

21.8 

23.3 

19 

3.4 

8.97 

- 13.7 

— 16 • 6 

- 18 . 1 

0.3 

15.72 

14.68 

11.45 

8.39 

5.85 

“ 1.7 

- 6.45 

- 8.2 

- 8.93 

- 7.0 

0.4 

8.97 

7.01 

m 

3 . 17 

1.95 

1.5 

B 

- 3.9 


- 3.8 

0.5 


m 

2 . 14 

2 

0.64 

0.78 

- 1.44 

- 1.5 

1.31 


0 - h 

2.47 

1.5 

0.8 1 

0.44 

0 

0.44 

- 1.2 

- 0.76 



O 

1 . 13 

0.6 

0 

0 

0 

0 





0.8 

0 

0 

0 

0 

0 



, 

1 



0.9 

0 

0 

0 
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n 2 £ 

ERRORS DUE TO LINEARISATION FOR S =0 ID 



0 

2 

0 

4 

0 - 6 

0.8 

1 

2 

4 

6 

8 

10 

0 1 



60 

67 

65 

' 63 

11 

40 61 

10 o: 

.-0 

48 

-7 8 

-12 6 

0.2 

g 




m 

24 33 

20 

69 

7 38 


-8 

14 

-10 2 

-11 2 

0 3 

g 





9.34 

7 

35 

0 64 

-3 5 

-4 

33 

-4 68 

-4 6 

0 4 

8 

95 

7 

ig 

5.41 

3 79 

2 

59 

-0 33 

gg 

B 


-0 91 


0.5 

4 

87 

3 

44 

2.42 

1 8 

1 

0 

B 

0 36 




0.6 

2 

35 

g 


1 . 02 

0.58 

0 

0 114 





0 7 

B 


0 . 

6 

gg 

0 18 

0 

16 






0.8 

0 

32 

-0 

16 

5.83 








0.9 












71 



























































TABLE 423 

ERRORS DUE TO LINEARISATION FOR S =0 iS 



0 2 

0 4 

0 

6 

0 . 8 

1 

2 

4 

6 

8 

10 

0.1 

40 6 

58 4 

64 

93 

65 . 53 

62 78 

IBB 


ID 

0.16 

-3 86 

0 2 

25 8 

30 

28 

18 

ID 


9 96 

0 17 

m 

BE 

- 5.33 

0 3 

SB 

ID 

M 

EB 

9 .99 

7 93 

2.2 

BQ 

-1 8 



0 4 

8 86 

7 39 

5 

D 


2 .95 

0 39 

-0 45 




ID 

4 65 

D3I 

2 

D 

SB 

m 

0 01 





0 6 


SB 

0 . 

89 

0 62 

0.27 






0 7 

1 06 

0 58 

1 

28 

2 , 35 







0.8 











0.9 
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TABLE 4.2.4 

ERRORS DUE TO LINEARISATION FOR S =0 20 


d'^ 

't,,. 

0 2 

0 4 

0.6 

0 8 

■■ 

2 

4 

6 

8 

10 


1 

40 

56 9 

63 3 

64 . 5 

63.09 

70 9 

22 6 

11 9 

IQQ 

5 88 


2 

25 5 

30 

28 5 

10.23 

22 

10.96 


-0 68 

QHj 

0 16 

0 

3 


14.8 

ED 

10.02 

8.1 

B9 

-0 3 




0 

4 

8,7 

7 .23 

5 44 

4.03 

29 

0 4 





0 

5 

4 6 

3 27 

2 25 

m 

0 93 






0 

6 

2 22 

1 34 

0.72 

0 3 

0 






0 

.7 

0.8 

0.22 

0 








0 

8 
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TABLE 425 

ERRORS DUE TO LINEARISATION FOR S =0 25 


D ^ 

0 2 

0 4 

0 6 

0 8 

1 

2 

4 

6 

8 

10 

0.1 

36.5 

51 3 

56 . 9 

58 6 


50 

30 

17 7 



0 2 

23.3 

26 8 

m 

24 5 

22 58 

BB 

IB 




0.3 




10 4 

8 8 

3 6 

0 16 




0.4 

8-08 

6.9 

5.33 

nn 

BB 

0 IS 





0.5 

4 28 

2.96 

1 . 94 

-0 3 

0 3 






0 6 


0.85 

0 12 








0 7 

0 21 










O.B 


































































4 3 EXACT SOLUTION FOR UNSTEADY STATE FLOW BETWEEN TWO 


PARALLEL DRAINS ON A HORIZONTAL SEMI -PERVIOUS BARRIER 
The height of drains above the impervious bed have been used 
effectively for linearisation for all times but when the drains 

are on the barrier this la not pos'^lble Realising 

this Glover (1954) found an approximate solution of governing 

partial differential equation without lineriaation The solution 

involved an integral to be evaluated numerically In the present 

study an exact analytical solution will be found 

431 NON-DIMENSIONALISATION 

The mathematical problem Is given by 


f) 


ah 

at 


= -s-C k b -s— ) 
dx^ s ax 


. ( 4.3 1 ) 


with 


h(x,0) = f(x) 


ax ( 0 , t ) 


= 0 


hCl^/2,t) = 0 

The co-ordinate system is shown in fig (4 3.1) 
Using the non dimensional variables, 


X = x/L 


H 

T 


= h/h 
' m 


K 


8 


t h 


m 




1 


2 

X 


(4.3.1 a) 
(4.3.1 b) 

(4.3.1 e) 


(4.3.2) 
(4 3.3) 

(4.3.4) 


eqn. (4.3.1) becomes, 



15 - £_ Th 1 

dT dH dX J 


wl th 


H(X,0) = F(X) 


i9U 


dX 


(O.T) 


0 


H(0 5,T) = 0 

where h la the aiaxiaum height o£ water table 
m ® 

the drains (1,6 at X=0)at T=0 
432 SEPERATION OF VARIABLES 


Letting, 


H(X,T) 


We get, 


1 d 

~T~ m 



] 


= F(X) G(T) 



= P 


where p la a constant Thus the following two 
equations are needed to be solved 


F f" + CF')* = p F 


with 

H3r (0) “ ° 


and 


F(0 5) = 0 


dG 

ST 


P 


G 


2 


With 


(4 3.5) 

(435a) 

(435b) 

(435c) 
midway between 

(4 3 6 ) 

(4 3 7) 
differ ent lal 

(4 3 8) 

(4 3.8 a) 

(438b) 

(4 3 9 ) 
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G(0) 


1 


(4 3 9 at) 


It should be noted that only h needs to be specified initially 

m 

as FCX) is found from the solution This is because 6qn.C4 3 6) 
indicates that profile at any time can be obtained from that of 
other time by multiplying it with a constant 
433 SOLUTION 

Eqn (4 3 8) can be written as 

FF* ^(Ff’) = pF® C4 3 10) 

Integrating, 

(FF )* != pF^ + a (4.311) 


and using eqn (4 3.8a) 



Sine F=1 at X=0 eqn (4.3 11) becomes 


(4 3 12) 


ff’ 



(4.3 13) 


Since 0<h<hm , 0<F<1 
expression under the 
Letting , 


and hence p must be negetive to get the 
square root sign a positive quantity 


p = - q 


(4.3.14) 


eqn (4.3.13) becomes 

- y/ q (1-F)(F* + F+1) 

Integrating , 


(4 3.15) 


F 


J 


F 


-OB -OB 

(1-F) (F®+F+l) dF 




<1 


X + b 


(4.3.16) 
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- J dF * j 


Ufllne, 


/y^ X ♦ b 


cos </> 


A - i + F 
A + 1 - F 


, A= -yTa 


C4 3 17) 


C4 3 18) 


and K* = sin* a = =0 9330127 (4.3.19) 
16 04=7 5® (4320) 
and noting that 


d 

dijfr 


6in4j6 y/ 1-K“sin“^ 


1+ cos^ 


1-COS^ 1 

1+cosiji yi-ksln*^ 


(4 3 21) 
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iollowlijg aolutioij o£ 6iitj (4 3 a) le obtained 


^ ^ -i i- aln ^ y 1 -K*sin*^ 

2 ( i) ECo»,<^) + C - X3) )F(«,<^)-2 (^r 

(*3)' 1 + coe 

2 

-3 q + b (4 3 22) 

when E and F are elliptical integral given by 

4* 

“ J / l-k'aln-* 

i' 

r -0 5 

F(^,a)= J (l-k®ein®^) 

0 

and b ia a constant of integration. Using X=0,F=l,<ji = 0, 
we get, 

b = 0 (4.3.25) 

Using X = 0.5, F = 0 ,4r = 74.457732 and tables provided by 
Abramowitz and Stegun (1944) to evaluate E (75,74 457732) 

and F (75, 74. 457732), we get, 

q = 4 4454231 (4,3.26) 

Therefore the solution for F(X) becomes 

2.632 E(oi,(^)- 0.556 F(oi,^)- 2 632sin«^ (l-k*aln*)’^ “(l + cos^)“‘ 

=1,7215156 X (4.3.27) 

Solution of (4 3 9) la given by 

1 

G = (4.3.28) 

q T + C 

Uith (4 3 9a) we get 

C = 1 (4.3.29) 

and hence G is given by. 


(4 3.23) 


(4 3.24) 
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G 


1 


(4 3 30) 


q"T + 1 

Hence the solution of eqn (4 3 1) is given by eqn (4 3 6) with 
eqn (4 3 27) and eqn (4 3.30). 

434 NON - DIMENTIONAL PROFILE 

The non-dimensional profile has been calculated by eqn (4 3 27) 
and shown in table (4.3.1) Profile at any time can be obtained 
by multiplying the h / values by a constant given by 

eqn (4 3.30) Profile ao obtained at various time steps has been 
shown in - 4 3 1 

Comparison with Glover’s solution, 

H = F(X) G(Y) (4 3 31) 

where F is given by 


X 


w 

r w. dw 



(4,3 32) 


and 


G 


1 

4.5 T + 1 


(4.3 33) 


shows that instead of 1 7215 and 4 4454, he has got y 3 and 
4.5 .Also this solution involves an integral to be evaluated 
evaluated numeriacally Note that in this solution origin is 
taken at left drain In his solution, Glover used initially 
horizontal water table, though the mathematical problem taken 
was same aa in eqn. (4. 3 1) From the curves provided by him, it 
can be seen that at T = 0.014 the profile la almllar aa given 
in table (4.3.1). 
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h /h 



Fig. 4.3.1 Profile at various times for flow 
between two drains on the barrier. 
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TABLE 4.3.1 


NON- dimensional WATER TABLE PROFILE 




X 

H 


H 

0 

1 



C.05 

0.9997 

0.30 

0.822 

'MO 

0.9981 

0.35 

0.6924 

0.1 5 

0.9989 

0.1 0 

0.5734 

0.20 

0.9082 

0.45 

0.4142 

0.25 

0.8524 

0.50 

0 


•o * M 


92 


4.4 LINEARISED SOLUTION FOR TRANSIENT SATURATED FLOW OVER 
SLOPING SEMI -PERVIOUS BARRIER 


The transient flow between two parallel channels constitutes 
one of the most important problems of gt ound watet flow Exact 
solution is non-existent Linearised solutions for different 
ca^es were found by Uerner (1957 ) , Glovet ( 1953 ) , Moody atid Tapp 
and Dumn (1960), Singh and Jacob ( 19 74 ) ,Sewram and Chuhan (1967) 
for flows over fully impervious beds with Dirichlet type of 
boundary conditions Lat inopoulas ( 19 86 ) considered third type of 
boundary condition at one end, the boundary condition at other 
end being either no flux Neumann type (impermeable boundary) or 
Dirichlet type for flow over horizontal impervious bed witli 
water table being horizontal initially Flow over horizontal 
semi -pervious bed was considered by Uessseling and Uesseling 
(1986) and they found a solution for flow between two drains 
when a steady rainfall occurs from top. In the present study a 
general solution will be found which is capable of handling all 
the three types of boundary conditions at either end for flow 
over a sloping semi-pervious bed 

4 4.1 NON-DIMENSIONALISATION 

The partial differential eqatlon governing unsteady flow over 
sloping semi-pervious bed is 


^h 


17 


St 


a 


k h 



dh 

( + I) 


+ + w 

c 


(4.4 1) 
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Uith initial and boundary condltlona aa 


h(x,0>=f(x) 


(4 4 1 a> 


-k 


s 


b 

1 


dh 




bkl+k h = kh atx = 0 

is i i 1 


C4 4.1 b) 


rfh 

k b +bkl+k h=kh atx=l (4.41c) 

sa. as a aa x' 

Sx 


Using Nondimens ional paramet ares , 


H = 


(4.4.2) 

X = 


. (4.4 3) 

U = 

w/k^ 

. . (4.4.4) 

tl 

k /k 

(4,4.5) 

II 

01 


. (4.4 6) 

T = 

kst/ni^ 

- .(447) 

c = 

o 

1* Ak.c) 

. (4.4.8) 

B = 

i 

b /I,. 

(4 4.9) 

B = 

a 

b /I 

ft X 

. (4.4.10) 

H = 
o 

Vx 

. (4.4.11) 

H = 

i 

V'x 

(4 4.12) 

H = 

a 


. (4.4.13) 


eqn (4.4.1) becoiaea 

aVL/dT = d/dX(HdH/dX) + I <>H/dX - C^H + + U (4.4 14) 
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With initial and boundary conditions as 

H(X,0) = F(X) (4 4 14 a) 

-E Sn/SK +KH=KH+BI (4414b) 

1 i 4 i 1 ' 

B <9H/dX + KH = KH-BI (4 4 14 c) 

2 2 2 2 2 

442 LINEARISATION 

As exact solution of eqji (4 4 14) la not poaeible the eqatloii 

la linearised following Glover (1953) as 

dH/dT = D H* * + I H*- CH4CH+U (44 15) 

o o o ^ 

with the Initial and boundary conditions i etnainitig unchanged as 

H(X,0) = F(X) (4 4 15 a) 

-B +KH=KH+BI 

1 i 114 

B ^H/dX +KH=KH “BI 

2 2 2 2 2 

443 TRANSFORMATION 
Terms In linear eqatlon containing dH/i9X and H can be 

removed through following transformation 


(4 4 15 b) 
(4 4 15 c) 


-(^iX+yT) T 

H = U6 + J'UdT (4.4 16) 


so that we obtain 


T (p X+w T) 

d\i/dT = D U" + ( ~ dT ) e .... (4 .4.17) 


with the transformed initial and boundary condition as 
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U(X,0) = F(X) 


SU 


-B 


(4 4 17 a) 


vT 


i dT~ ^ ^ S + I^l) = [ Bi. I - /U dT ] e 

(4 4 17 b) 




[ K H - B I 
2 2 2 


K dT ] e 

2 0 


i^T 


(4 4 17 c) 


where 


^ = I/2D (4 4 Ifl) 

w = D + (4.4 19) 

The nonlinear term in eqn.(4 4 17) can be removed tlirough a 
further t rana f ormat ion 

C/liX+wT) 


U = V + V e 


(4.4.20) 


ao that we get 


dV/^T = D V" (4 4 21) 

where y is given by the following diffrential eqation 

T 

= Dv'**+IV*-Cv+CH-C J“UdT (4 4 22 ) 

The solution for y/ depends on the form of^H (X,T) and UCT) 
and can be found by assuming y(X,T) to have a similar form as 
H^(X,T) in X 
Thus for 

H = P X + Q . (4 4 23) 

o 

asssuming v^(X,T) as 

y; = T^(T) X + T^CT) (4 4 24) 

and putting in the above eqation for v' we get two differential 
equations as 

dT /dT =-CT+CP (44 25) 

* o i o 

T 

dT /dT = I T - C + CO - XU dT- (4.4 26) 

2 ' i O ft O O O 


B6 



as 


Solving eqn (4 4 2S) we get T 

i 

T = P 

1 

Using this eqn (4 4 26) becomes 

T 

ST /ST + C T = P I + C Q - C /U dT 

® 0 2 O O O 

T =6 KPI + C Q-C /u dT)e^o^ dT + C ' ^ 

2 V 0€»0 


(4 4 27) 


(4.4.28) 


(4 4 29) 


Selectine arbitrary conatant C as zero and approximating 


recharge as 


U(T) = 2 T + bj 


j=i.2.a 


T < T < T 
^ ^ ^ j + 1 


J=l,2,3. n 


T = 0 

1. 


T = T . 
j+i end 

the expression £or T becomes 

a * PI+C Q- & 

1 2 €3 

T = -Jt +o-T + 

2 u U 

. a-C b* 

where j o j 


(4 4 30) 


(4 4.31) 
. . (4 4.32) 


(4 4 33) 


(4 4 34) 


Thus finally we obtain, 


SV/ST = D S^W/SX^ 

with initial and boundary conditions given by 
V(X,0) = [F(X)-v'CX,0)]e*^'^^^ 


(4.4 35) 


(4.4.35 a) 


■^- 2 ^ ^ 


B )V = ,Jk I -K, ,/U dT -K v;( 0 , T ) + B, 6^“^ 


(4.4.35 b) 




fk H 
]_ a 2 


I -K /U dT “K Cl 
2 2 0 2 


,T) Je 

(4 4 35 c) 
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with the tr ana forming function 


V = (H-XUdl-v^) giJX+vT 
The backtrana format ion ia given by 


H = V e + yu dT + w 

U 


From eqn (4 4 24) 


dv-' 

ix 


= p 


Thua finally we want the aolution of , 

^ ^2v 




d X- 


with initial and boundary condltiona given by 


V(X,0) = t F(X)-v'(X,0)]e^ 


-B 




+ E V = G (T) 
dX 


c 

(4 


B + E V = G (T) 
2 — ■ 2 2 
dX 


(4 


where 

E = K + p B 

i * 1 

E2 = K2- B2 

GCX)= ( F(X)~%y(X,0)]e^^ 

’ vT 

G = (K H + B^I - K dT-K^ V'(0 ,T) + B^ P)© 

* * ^ ' vT+u 

G = (K H - B I - K /U dT-K y(l,T)-B P)e 

The aolution of above problem ia given by , 


(4 4 36) 

(4 4 37) 


(4 4 38) 

(4 4 39) 

I 4 39 a) 
4 39 b) 

.4.39 c) 


C4.4.40) 

(4.4 41) 
(4.4 42) 

(4 4 43) 
4 4.44) 



V 


(4 4 46 ) 


where 


where 


N=i 


B Y 

'NX e N ( ^ “n 


is the root of 


tan a 


c*-, ( B E +B E ) 
N 12 2 1 


N 


(4 4 46 ) 


ft 


B B - EE 

1 2 N 12 

2 b“ £*„ + t 
2 N 2 


N 


CB®01^, + E^)[(B a®+Ef) + B E ]+B E (Bo®, +E ) 
'iN I' l N 1 22 11 ' 2 N 2 


y^v ® B o(„ Cos o„X + E Sin ot„X 
NX 1 N N 1 N 


(4 4 47 ) 
(4 4 48 ) 


= J G(X) dX ( 4 . 4 . 49 ) 

o 

’'NT “ / 
o 


i 

C = J G(X) dX 
o 


= 0 


If 


E =E =0 
1 2 


Otherwise 


( 4.4 51 ) 

(4 4 52 ) 


|. ““ ““Vj - “n ■ V 

n,=-5^Ca„ sin o,j + cos 0.^-1 

N «N 


(4 4 52 ) 

(Oj^cos Ojj+ sin Ojj) 
(4 4 54 ) 


M 

9 

M 


4 


L M - L n 
11 9 2 

L M + L M 

fi 1 9 2 


(4 4 55 ) 
. (4 4 56 ) 


L = 


E (B +E )/ L 

12 a 4 


L = 


B E yL 

12 4 


L = 


E E /L 

12 4 
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( 4 . 4 . 57 ) 

. (4 4 . 58 ) 
(4 4 . 59 ) 


L = B 


Fl+ B, E^+ E, 


(4 4 60) 


To evaluate the initial condition integral the initial 

condition ia approximated by a set o£ parabolas each fitting 
three consecutive points of F(X) as 


ECX)=2 X® + b^X + c^ for X^:S X ^ X ^^2 ^ ^ 


i=i-a,E 


X = 0 

i 

X = 1 

m+2 


(4 4 62) 
(4 4 63) 


E ( [a^X 4C b^-I,)X + c J + E I)-[2a X + b -T ] 

i=* ,fi ,a 

I ^i+a 

rj +E^ J^I ) + [2a^B^«jj XXJ + 2a^E^ XXI ])l 

' X 

^i 

(4 4 64) 

where 


I = /Sin o(„X e^ dX 
N 


(4 4 65) 


J = /Cos OjjX ^piX 


(4 4 66) 


I = - 


2^2 
«N + 


( otj^ cos otjjX - sin <^jjX) ( 4 4 67 ) 


j ( ^ cos oij^X +«jjein otj^X) (4 4.68) 


JI = - 




J + 




(4.4 69) 
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Jj = - 


01 


N 


"n * 




oi 


N 


^ 2 ^ 2 


/J J =- 


'N 


2^2 
"n ^ ^ 


SJ + 


JI + 


2 

+ 




2 

fJ 




2^2 
«N + ^ 


/'I 


(4 4 70) 


C 4 4 71 >l 


(4 4 73l) 


To evaluate the boundary condition integral, the transient 
boundary condition is approximated as 


H^CT) = 2 C Cj T + d^) 
H,(I) = 2 « ‘‘j I * fj) 


T <T<T 
J ^ J + 1 


(4 4 73d.J 
(4 4.73 b) 


j = 1 , 2 , 3 . n 

Then, the boundary condition integral cati be writteti as 

n 


"ni'2 


(4 4 74) 


.l=i 


where 


kT 


I » [A T+B -« A ] 

4^1 i. 4 


(4.4.75) 


M 


I . [A T+B, Ao ]- 

ft '■ Z 2 S 


«T 


(4.4.76) 


M 
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A = 

B 

B = 


K (-b +c -cr ) 

* J J 

K (-b +c - 0 -) 
a J J 

K (d - 

* J 


= K (f - 
a J 


PI+C Q-<y 

Cl 


PI + C 0- «y 
o 


and 


= D a 


N 




Thus boundary Integral can be evaluated 


) + B^(I + P) 
)-B (I+P) 


(4 4 77) 
C4 4 78) 
(4 4 79) 

(4 4 80) 

(4 4 81) 


445 APPLICATION TO A DRAINAGE PROBLEM 

The solution obtained above has been applied to a problein 
involving flow over a sloping semi -pervious bed. The data and 
analytical results has been shown in /ftf-4 4 1 
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Recharge from top 
w= 0 0A2 e”* m /d 

1 i i i i I 



Fig. 4 4.1 Flow between parallel drains over sloping 
semi -pervious bed 


<13 





4.5 LINEARISED SOLUTION FOR UNSATURATED FLOW ABOVE WATER TABLE 


To aasea the recharge to ground water table due to rainfall 
at ground surface requires the solution of Richards’ equation 
for unsaturated flow Moat of the solutions found till date are 
for a semi inf inlte soil column with uniform initial saturation 
Some of them consider infiltration when the surface of soil is 
saturated whereas others assume a constant late of Infllttatlon 
at ground surface 


451 GOVERNING DIFFERENTIAL EQUATION 

The partial differential equation governing the transient 
unsaturated flow is 



^ ( 

. . s& . 




— ) - 

St 

Sz 

dz 

Sz 


( e^e + e^ ) 


(4 5.1) 


with initial condition 

d(z,0) = f(z) 

and boundary conditon at water table, z = 1 as 

Z 

= a, 

and boundary condition at ground surface, z = 0 as 


(451a) 

(451b) 


t < t 


Si 


- A ^ + k = r 


t < t < t . 

Si di 


e = 6 


‘d.< <''b2 


- 6?S. . k = 

Az 


.(4.5 1 c) 

(4.5 1 d) 

(4 5 1 e) 
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where 


& volumetric moisture cont ent , non-dlmens ional ^ 

^8 saturation moisture content , non-dimensional , 

z distance measured from ground surface vertically 
downwards , 

1^ distance from ground surface to water table, 

6 dlffusivity of soil mass.L®!”*, 

unaaturated permeability of soil mass.LT"*’, 

T actual transpiration rate per unit depth of root 

a 

2one,T~*^ , 

r rainfall rat6,LT~*j 

e potential evaporation rate.LT"*; 

t time since start of rainfall,!, 

t^^ time of beglnlng of first rainfall,!, 

t time to saturation of soil surface since the beginning of 
rainfall ,!, 

t^^ time when drying of soil surface starts since the beginning 
of rainf al 1 , T , 

^bz time to start of second rainfall since the beginning of 
firat rainfall,!, 

Using the non-dimensional parameters 




e 





(4 5 2 > 


Z = s / 1^ 


( 4.5 3 ) 
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k t 

B 

1 c& - & ) 

z a r 

6 ( s - r ) 

k 1 
a z 


(4 5 4 ) 


(4 5 5) 


C4 5 6) 


R = £- 


(4 5 7) 


E = _ e C de - dr ) 

^ k ^ 
s 


(4 5.8) 


where 


_ c ) 

8 


residual moisture cont ent , non-dlmens ional , 
k saturated permeability of soil masa.LT'* 


the mathematical problem becomes 


(4 5 9 ) 


^ = £_( A 

dT SZ dZ 


- + E^) 


(4.5 10) 


with initial condition as 


,;6(Z,0) = F(Z)_ 


(4.5 10 a) 


d - d 
s r 


and the boundary condition at water table, Z *= 1 as 


4^(1, T) = 1 


(4 5.10 b) 


and boundary condition at ground surface, Z = 0 as 
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T. < T < T 

bl Bi 


- A ^ + K = R 


(4 5 10 c) 


T < T < 

di 


0(0, T) = 1 


(4 5 10 d) 


T 

Si 


< T < T 


di 


- ^ t K 

az 


E 


.(45 10 e) 


4 B 2 TRANSFORnATION 


An exact solution of eqn (4 5 10) subjected to eqn (4 5 10a) 

eqn (4.5.10 b) and one of the eqn (4.5.10 c,d,6) is difficult 

because of its non-linearity and hence a linearisation approach 

becomes inevitable Such approaches have been used earl 1 et , which 

consist of replacing A by A and K ^ by K ^ but were found 

av az ^ Ql-v az 

to be of limited use 

Eqn (4 5 10) can be written as 


wher 6 


Sd) ” * ' 2 ' ' 

^ = A0 + A0 -K0 - - E- 

ar 1 ^ 

(4.5 11) 


, a^ 

^ = (4.5 12) 

az 

* " TT^ 
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(4.5 13) 



t 

A 

t 

K 


<JA 

dK 


C4 5 14) 

(4 5 15) 


Replacing A by A 

n A V 

t t ^ 

A <f? zero Instead 


as used by earlier inv es t iftators makes term 
we can use following transformation 

4. - F(U) ‘ “) 


ao that the transformed equation becomes 


<»U 

dT 


A U - K U - ^ ( + E, ) 

S4, ^ ^ 


with F(U) given by 

H f t >} 

AF + AF =0 


(4 5 17) 


, (4.5.18) 


Using following K~<it> and A— 0 relationships 

K = . (4 5 19) 

and 

A = (4 5 20) 

Eqn (4 5 18) can be solved and a particular selection of 
integration constants appearing in the solution yields 

^ = --^1 (4 5 21 ) 

U 2 


SO that the 


From eqn (4 


U = ^ 




eqn. (4 5 17) now becomes 
dU 


9 § 9 9 


^- = A U - K U -CA^ + DE^U -(A2 + 1)E^ 


1 22 ) 


dU . 

_ = (A^. 1) 4 > a 


. (4 5 22) 

(4.5.23) 

, (4.5 24) 
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Writing A^K and K in terms of U 


A - A U 


(fl 1 25) 


K = K U 
1 


K x'CA +1) 

ft 2 


k'= K K 
1 e 


C4 1 26) 
C4.5 27) 


and therefore the eqn (4 5 17) becomes , 


JT * d 2*^ 


- K K 


i. s 


dU 


-C4k + 1)EU-(A +1)E U 


with initial condition ae 


UCZ,0) = [ F(Z)3^2‘^^ 


(4 5 28) 

(4 5 29 a) 


and the boundary condition at water table, Z = 1 as 

UC1,T) =1 . (4 5 29 b) 

and boundary condition at ground surface, Z = 0 as 


< T < T 

hi. SI 


t dU 


(A^+1) dz 


♦ K U ' R 

4 


(4.5 29 c) 


T < T <T . 

Si di 


U = 1 


(4 5.29 d) 


T < T < T. 
Si b2 


+KU® =-E 






(4.5.29 6) 
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The advantages of tanafotmatlon given by eqn (4 5 16) are 
now evident The term containing A' has been taken into account 
Also now exact and linearised solutions of eqn (4 5.29) will 
differ leas because powers of ’U’ requiring linearisation are 
genetally less than 1 Uhen defining back the profile in terms 
of </> , U will have to be raised by power of t— T i < 1 always 

A T 1 
B 

and hence the errors will further go down approximately by a 

by a factor of j- 

A +1 
a 

4 S 3 LINEARISATION 

The transformed unaaturated flow equation can be 
so that the linearised equation becomes 

~=L U"-L U'-L U-L 

dT ‘ * fl 4 

with initial condition as 

A +1 

U(Z.O) ={FCZ)} “ 

and the boundary condition at water table, Z = 1 as 

U(1,T) = 1 

and boundary condition at ground surface, Z = 0 as 


linearised 

(4 5 29) 

(4 5 29 a) 

(4 5 29 b) 


T 


bi 


< T < T 

81 



L U 

d 


R 


(4 5 29 c) 


T < T < T . 
Si di 


T 


di 


< T < T 


be 


where 


U = 1 


•L ^ t L U 
- dZ ^ 


- E 


2 


L = A 

i i 


(4.5 29 d) 

(4 5.29 e) 


. (4.5.30) 
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K -1 
2 

A +1 

L = K K U ® (4 5 31 ) 

a 1 ft 


L = (A +1) 

a ft i 


L = (A +1)E U 

4 8 2 


A +1 
2 


L = 


L = 


A +1 

ft 


i. 


K,-CA^-l) 

K U 


C4 & 32) 

C4 5 33) 

(4 5 34) 

(4.5.35) 


Once the solution In terms of U is found 4 > can be obtained 
by backtrana formation given by eqn.(4 5 21) and recharge to 
water table will be given by 


U 


K 


A +1 

B 


dU 

sz 


+ u 


A +1 
e 


at 2=1 


(4 5 36) 


454 SOLUTION 


To obtain a aoluiton of eqn (4 5 29) following subatituition 
is used 


-pZ-wT 

V = (1 - U + \y) e 

where (j , w and y selected such that , 

2L 

i 

V = “ L 


. (4.5 37) 


(4.5.38) 


(4.5 39) 


and 


y. = (- + L^)/L^ 


. (4.5.40) 
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U In terms of V will be given by 


U = 1 - V e 

The eqn (4 5 29) now becomes 


/jZ + wT 


(4 5 41) 


i>V 

— = L V" 

sr * 


(4 5 42) 


with initial condition given by 


A +1 

V(Z,0) = [1 -(F(z)) ® 1 e 


(4,5 42 a) 


and the boundary condition at water table, Z = 1 as 


V(1,T) = v> 6 




(4 3 42 b) 


and boundary condition at ground surface ^ Z = 0 aet 


T. ^ < T < T 

bi Si 


L ^ +( hL - L )V = [R + L ^ - L V - L ]e 

“ Jin B <S B 3Z <f ^S 


.(4.3 42 c) 


T < T < T . 

SI di 


V(0,T) = y; e 


(4 3.42 d) 


T . < T < T. 
di bi 


L5 — +(^L - L )V = [-E + L - Lw 

B iS ^ €F ^ ^ 


L ]e 

(4 3 42 e) 


Note that In this case 


(4.3.43) 


Thus finally we obtain the equation and initial and boundary 


tonJtt 


mt o.i 




. , . . C4 i,44 ) 


With Itiitial eotiditioh Ati 


V(Z,0) = G(Z) 


(4 3 44 a) 
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and boundary condition at water table, Z = 1 as 


V(1,T) = G (T) 


(4 3 44 b) 


and boundary condition at ground surface, Z = 0 as 


T. < T < T 
bi sr 


+ L V = G (T) 


(4 3 44 c) 


T < T < T 

Si dl 


V(0,T) = G ^^(Tl 


( 4 3 44 d 1 


^d.< T < Tb. 


av 

‘■= -5Z- * S'' = G,,(T) 


(4 3.44 el 


where 




(4 3 45) 


A +1 


G(Z) = [1 -(F(Z)) + v' ] e 


G^(T) = v e 


-vT-fj 


(4 5 46) 
(4.5.47) 


G^/T) = [R + - L^l« 


-wT 


(4 5 48) 


G^^ (T) = V e 


-vT 


(4.5 49) 


dy/ 

G^3(T) = [-E + ^ - L^]e 


-wT 


(4 5 SO) 


Solution of eqn (4 3 44) during the three phases is obtained 
as follows 

PHASE-I T. < T < T 
bl at 

Solution of eqn (4.5 44) subjected to eqn (4 5.44a), 
eqn (4.5.44 b) and eqn (4 5.44 c) is given by 
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oa 


V = 


I 


-la T 


(y y e 


^ I- r o n 

^ ^ I 2 I 

+ L a T I 
N i H KT I 


where is the root of 


and 


L tan oi - L ot 

N B N 


(4 5 Bl) 


(4 5 52) 


" L L 

B 7 


(4 5 53) 


= - L « coa ot Z + L sin a 2 

MB B H M 7 H 


MT 


M = 


I 

a 

T . fit 

r L a X - 

f-" i' 


G(Z) r dz 

MB 


G (\) + 
a ii 


M, G^(N5] 


Ck 




L sin oi - — ( coa oi 


w- ] 


C4 5.54) 
(4 5 55) 


(4.5 56) 
(4 5 57) 


M 


= [- 


L L , 

B - 7 1 

~ ( sin + cos 0 (^,-l) “ — ^ (oi^^cos + sin 


ot N 




N 


W 


CA 


2 N >4 


”a = (L - L ) 

7 B 




[ H - W ] 

i R 


”4 = (L - L ) t -VM,+ L,I1,] 
'7 k' 


(4.5.58) 
(4 5 59) 

(4.5.60) 


The two Inteerala 4?^ and are respectively called as 

initial and boundary condition integral They can be evaluated 
as follows 

A +1 

Approximating the expression [1 -(F(z)) + V ] 
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[1 -CF(z)) 


(4 5 61) 




■^ V']=y afz + bZ+c 


^ = 0 
^ n +2 ~ ^ 


2 ^ ^ 1.2 

1 = 1,3,5, m 

(4.5 62) 

(4.5.63) 


so that G(Z) bacomea 


IVi 


G(Z) = [ J a^z . b ^2 * cl 

^ 1 = ITb , i 5 J 


+ c I e (45 64) 


Z ^ 2 :S 2 
1 1 + 2 


Uaing thia the Initial condition integral becooiea 
Z 


=N = 1 

L~i,a. 


1+2 


J (a 2"^+ b.Z + c,)e ^^(-L oi coa a Z + L sin a Z)dZ 

•^1 i i BW NV M 


B Z 


i=?a,B [ ' "m ] 


(4.5 65) 
(4 5.66) 


where 


I = J Ca|Z^+ b^Z + c^)e"^^ ain a.Z 


N 


(4.5 67) 


Ca^Z^ b^Z + c^) - (2a^Z+bp Jl^ + ZaJJl, 


(4 5 68) 


J = 


= J ^^i 


Z % 


b,Z 


C ^)6 


-{JZ 


coa a. Z 

N 


(4.5 69) 


= (a^Z^ bjZ * op 3 ^ - (ZajZtbp ■* 2ajp^ 

(4.5 70) 
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(4 5 71 ) 


■c = I 


- I 


Ih 


p, 


;p, 


jp= = 


e sin 

a Z 

M 

dZ 



L 

cc»a 

ot 2 + /- 

2 , 2 
d 4J 

KI ^ 

r» 

-uZ 

6 COS 

« z 

H 

dZ 


e-^2 

- [« 

sin 

7 - 

2 2 
d + U 

N ^ 

L « 

Ol ^ 

Ol 

w 

.7 


t-i 

2^2 

Oi + /J 

li 

C 

2 

d 

L * 

+ M 

d 

T 


AJ 

2 ^ 2 

a + 

X 

C 

2 

d 

a. 2 

+ /J 

d 

N 

Pc 



2 . 2 
d + /J 

2 

d 


d 

fl 


AJ 


(4 5 72) 


- /j cos a 


(4 5 73) 

(4 S 74) 

C4 5 75) 


(4 5 76) 


Jl • . (45 77) 


c 2 . 2 *v c 

d + p a ft 

^ Ol + p 


Pc ^ 


To evaluate first I and J and their integrals 


/J"! and SSjc are evaluated from whiul> I ai»ci J Cfi^h ^'PWnd At 

<1 

Z, and 2 , Then from eqn (4.5 66) can be found 

i i z 
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To evaluate the boundary condition integral rainfall can be 
wr itien as 


RCT) 




a T+b 
J J 


(4.5 79) 


T <T<T 

J J + i 

J = l,2,3 ri 


T = 0 
1 


(4 5 80) 


■'nt. ' ■'s, S- 81 ) 

Using this approximation for rainfall tha boundary condition 
integral becomeB 

T . a 
La K 


MT 


r Li a - 


I 

= [ ". 'b. " " 4 I- ]. 


(4 5 83) 


where 


Ibi f ^ [R-L^(l+V)] d\ 


(4 5 84) 


ot: 


«y { a T + b -L (1+v)) - a 

ft I ' J J a'' 


J 


. (4 5 85) 



v 

«yT-/j 

A ~ 

bft 

& 

"• t3 

<y = 

L 

2 

Oi - V 


1 N 


(4 5 86) 
(4 5,87) 


Thus by evaluating I, and I, can be evaluated. 

04 02 NT 


PHASE- 1 1 T < T < T. 

81 <il 

Solution of eqn (4 5 44) subjected to eqn.(4 5 44 a), 

eqn (4.5 44 b) and eqn (4.5.44 d) ie given by 
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<30 


V = 2 


£ S 

7 " ■ 


SinCNnZ) [ £ + L Nn T ] 

M i NT 


N= 1 .S,3 


wher e 


(.4 . & ob; 


= J G(Z) sin (NnZ) dZ 


(4.5 89) 


and 


L N^n^T M 

NT J - ‘ [G^(T) - (-1)” G^(T)3dT 




(4.5 90) 

Again approximating G(Z) aa in eqn (4.5 64) £ is given by 




* = \ 


[I]. 


'i + 2 


(4 5.91) 


where I is given by eqn (4.5 66), Boundary condition Integral 


T is given by 

NT “ ' 


T = y - f-n” (L 

LN“n“-N 

i 

-1 + (-I)"" e"^ ] (4 5 92 ) 

455 VALIDITY OF RESULTS 

Linearisation of the governing partial differential equation 
introduces error in the predicted moisture content 

Consider a problem in which the soil surface is saturated and 
infiltration is taking place in the soil which is initially at 
uniform moisture content. As the water will inf i Itrat e , more and 
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more volume of soil maas will get aaturated Unaat ui: at ed flow 
will continue till the entire aoil maaa geta aatutat ed . For thia 
moisture content profile at various times has been calculated 
by analytical as well as numerical approach and the results are 

shown in 5 1) The data used for calculation for the soil 

considered are 

\ = 1 8 

P, = 75 
b 

d - e = 0 27 
r <iip 

& = 0 47 

s 

= 0 217in/d 
0^ = 0 37 

The data has been taken from Jensen and Hanks (1967) and 
Brooks and Corey (1966) for a silty loam From the /i-tf“(4 5 1) 
it can be seen that the moisture content upto wetting front can 
be predicted nicely beyond which it is overestimated Linearised 
solution overestimates moisture content at all times and this 


over estimation increases as one goes fatther from this point. 
The remarkable accuracy with which the linearised solution 
predicts the time to saturation for a particular point can be 
noticed from the graph The maximum error in the effective 


saturation (0-0^ ) was 53% whereas in the volumetric 



moisture content it is 14 54 % The recharge rate to the water 


table was also calculated by both the methods but errors due to 
linear iSAtiotj wet w found to be unacceptably large(500%) This is 
because nonuniform error in 0 will result in large errors in ^ 
and hence in recharge rates 
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0 



Fig. 4.5,2 Flow in unsaturated zone when 
surface is saturated. 


no 




CHAPTER B 


CONCLUSION 

B . i I NTRODUCTI ON 

In this study a fisneial 1-D ground water flow problem in 
rectanguler coordinates was considered 

Exact solution for steady state flow over a horizontal 
semi-pervious bed was found by using the solutiotis provided by 
Sikkema and Van Dam( 1 982 ) ,N1 euwaal and Ni Jhuia ( 1 9 7 9 ) As exact 
solution for the case when bed is sloping is difficult to find, 
a linearised solution was found for this case and effect of 
linearisation on accuracy of prediction of water table heights 
and location of ground water divide was studied 

For unsteady flow between two parallel drains provided on a 
horizontal impervious barrier, an exact analytical solution was 
found and profile of pheratlc surface was computed 

For general transient flow problem in finite domain exact 
solution is not available and hence a linearised solutions is 
found The solution is applicable for arbitary initial condition 
and transient boundry conditions of any of the three types 
viz Dirichlet,no £ lux(Neumann) and Cauchy, for flow over sloping 
semi-pervious beds Substitutions have been used to remove most 
of the terms of the llnrearlsed equation and the transformed 
eqatlon so obtained la a simple one for which solution can be 
easily found using the method of separation of variables Subs- 
tistutions make boundry conditions transient. The solution for 
transient boundry conditions has been obtained by applying 
Duhamel's theorem 
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B. 2 GENERAL CONCLUSION 

Ftoiu the above study followine conclualons can be made 

(1) The Unerased solution for flow between two parallel drains 

over horizontal a emi-pervloua bed compares well with Young’s 

solution with errors generally less than 5% for all d/h i 0 2. 

o 

This suggests that linearised analytical solutions can be used 
ef f ectlvely 

(2) The linearisation approach is valid as long as d/h^> 0 7 as 
errors will be leas than 5% . Instead if 10% errors are allowed 
the d/hm greater than or equal to 0 6 can be considered 

(3) The linearised solution for unsaturated flow predicts the 
moisture contetit fairly well till the wave front but thenafter 
it overestimates the moisture content 

(4) The lineal ised solution for unsaturated flow cannot be used 
for predicting recharge to water table 

6.3 LIMITATIONS OF THE STUDY 

The solutions found in the present study at e valid only for 
a homogeneous and isotropic non-def ormable porous media Darcy's 
law has been taken to be valid for both saturated as well as 
unsaturated flow The artesian head in the underlaying aquifer 
should not be affected appreciably due to leakage to unconflned 
aquif er 

5.4 SUGGESTIONS FOR FURTHER WORK 

(1) The solution based on linearisation approach was found 

valid for d /h> 0 7 Efforts may be made to find better results 
ro 

for d/h <07 also 
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(2) The solution found here can be used to study the effects of 
the clogged banka and transient boundry conditions on water 
table heights 

(3) The substiutiona used in this study should be tried for 
obtaining solution for 2-0 flow between rectangular boundaries. 

(3) Various forms of and p-<^ relations should be investi- 

gated to find the transforming function of the unaaturated flow 
equation so that a proper form tan be selected for accurate 
prediction of soil moisture content and water table recharge. 

(4) Finally the applicability of above approach should be 
be investigated to the multiple aquifer problems 
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